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Abstract 

We provide a simple analytic formula for the two-loop six-point ratio function of 
planar A/" = 4 super Yang- Mills theory. This result extends the analytic knowledge 
of multi-loop six-point amplitudes beyond those with maximal helicity violation. 
We make a natural ansatz for the symbols of the relevant functions appearing in the 
two-loop amplitude, and impose various consistency conditions, including symme- 
try, the absence of spurious poles, the correct collinear behaviour, and agreement 
with the operator product expansion for light-like (super) Wilson loops. This in- 
formation reduces the ansatz to a small number of relatively simple functions. In 
order to fix these parameters uniquely, we utilize an explicit representation of the 
amplitude in terms of loop integrals that can be evaluated analytically in various 
kinematic limits. The final compact analytic result is expressed in terms of classi- 
cal polylogarithms, whose arguments are rational functions of the dual conformal 
cross-ratios, plus precisely two functions that are not of this type. One of the 
functions, the loop integral Q!^'^\ also plays a key role in a new representation of 
the remainder function TZ^^ in the maximally helicity violating sector. Another 
interesting feature at two loops is the appearance of a new (parity odd) x (parity 
odd) sector of the amplitude, which is absent at one loop, and which is uniquely 
determined in a natural way in terms of the more familiar (parity even) x (par- 
ity even) part. The second non-polylogarithmic function, the loop integral VL^'^\ 
characterizes this sector. Both fi*^^) and fi*^^-' can be expressed as one- dimensional 
integrals over classical polylogarithms with rational arguments. 
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1 Introduction 



Much progress has been achieved recently in the analytic understanding of seemingly com- 
plicated scattering processes. In particular, attention has been focused on the planar sector, 
or large N limit, of maximally supersymmetric JV = 4 Yang- Mills theory. The scattering 
amplitudes in this sector of the theory obey many startling properties, which has led to the 
hope that the general scattering problem might be solvable, exactly in the coupling. 

One of the major simplifications that the planar A/" = 4 theory enjoys is dual conformal 
symmetry [H El |3l IH El El [7] , which dictates how colour-ordered amplitudes behave under 
conformal transformations of the dual (or region) variables defined via Pi = Xi — Xj+i. For 
the particular case of maximally helicity-violating (MHV) amplitudes, this symmetry is 
intimately connected to the relation between the amplitudes and Wilson loops evaluated on 
polygons with light-like edges, whose vertices are located at the Xj [6l El El HI El El UHl [11] . 
The MHV amplitudes are infrared divergent, just as the corresponding light- like Wilson 
loops are ultraviolet divergent. The Wilson- loop divergence has the consequence that a 
suitably-defined finite part transforms anomalously under the dual conformal symmetry [H 
E]. The Ward identity describing this behaviour actually fixes the form of the four-point 
and five-point amplitudes to all orders in the coupling, to that given by the BDS ansatz |12j . 
From six points onwards, the existence of dual-conformal invariant cross-ratios means that 
the problem of determining the MHV amplitude reduces to finding a function that depends 
only on the cross-ratios — the so-called 'remainder function', which corrects the BDS 
ansatz. 

Great advances have been made recently in understanding the form of the remainder 
function, which is non-trivial beginning at two loops. The need for a two-loop remainder 
function for Wilson loops was observed for a large number of points in ref. [7], and for six 
points in ref. [9]. The multi-Regge limit of the six-point scattering amplitude also implied a 
non-trivial remainder function [13]. At a few generic kinematic points, the Wilson loop [11] 
and amplitude [TU] remainder functions were found to agree numerically. The six-point Wil- 
son loop integrals entering the remainder function were computed analytically in terms of 
Goncharov polylogarithms p^lTS], and then simplified down to classical polylogarithms |T6] 
using the notion of the symbol of a pure function [T71 [181 (El HB] ■ The integrals contribut- 
ing to the six-point MHV scattering amplitude have also been evaluated analytically [20] 
in a certain kinematical regime using a mass regulator [2T], and the remainder function 
has been found to agree with the Wilson-loop expression of ref. [IS]. Very recently, the 
symbol for the three-loop six-point remainder function was determined up to two arbitrary 
parameters [22], by imposing a variety of constraints, in particular the operator product 
expansion (OPE) for Wilson loops developed in refs. [231 1211 [25] . 

For more than six points, numerical results for the remainder function have been ob- 
tained via Wilson loop integrals [26] [27]. Integral representations for the MHV amplitudes 
have been presented at seven points [28] and for an arbitrary number of points using mo- 
mentum twistors [in]. Recently, an expression for the symbol of the two- loop remainder 
function has been given for an arbitrary number of points [30], and the structure of the 
OPE for this case has been explored [31]. In special kinematics corresponding to scattering 
in two space-time dimensions, analytic results are available for a number of configurations 
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at two loops [32| l33t [M] and conjecturally even at three loops [35] . 

When one considers amplitudes beyond the MHV sector, there is another finite dual 
conformally invariant quantity that one can consider, namely the 'ratio function' V. This 
quantity is defined by factoring out the MHV superamplitude from the full superampli- 
tude I36j, 

A = AMYiV X V. (1.1) 

Infrared divergences are universal for all component amplitudes; hence the MHV factor 
contains all such divergences, leaving an infrared finite quantity V. One of the central con- 
jectures of ref. [36] is that V is also dual conformally invariant. There is strong supporting 
evidence for this conjecture in the form of direct analytic one-loop results [37 t [38 | [39lH0t 1^ 
and, in the six-point case, numerical evidence at two loops [32]. In this paper, we will con- 
struct the ratio function V analytically at two loops for six external legs. 

At tree level, the ratio function is given by a sum over dual superconformal 'i?- invariants' 
[36| H3| HU HSl HHl H?]. These quantities are invariant under a much larger (infinite- 
dimensional) Yangian symmetry, obtained by combining invariance under both the original 
and dual copies of superconformal symmetry [IS] . Beyond tree level one finds i?- invariants 
dressed by dual conformally invariant functions [36llll]|l2]. Since the i?- invariants individu- 
ally exhibit spurious poles, which cannot appear in the final amplitude, they cannot appear 
in an arbitrary way. The particular linear combination appearing in the tree amplitude is 
free of spurious poles. At loop level, the absence of spurious poles implies restrictions on 
the dual conformally invariant functions that dress them [19|. Additional restrictions on 
these same functions come from the known behaviour of the amplitude when two of the 
external particles become collinear. These constraints will be important in our construction 
of V at two loops. 

A consequence of the duality between MHV amplitudes and light-like Wilson loops is 
that the remainder function can be analysed by conformal field theory methods, such as the 
operator product expansion (OPE) [231 l2ll [25]. Various proposals have been put forward 
for extending the duality between amplitudes and Wilson loops beyond the MHV sector, 
either in terms of a supersymmetric version of the Wilson loop [501 [511 [30] , or in terms of 
correlation functions [521 [S3]- Although there may be various subtleties in realising such an 
object, compatible with the full A/" = 4 supersymmetry in a Lagrangian formulation [53] , 
one may instead justify the existence of such an object through the OPE. The framework 
for pursuing this approach was developed in ref. [55], and agreement was found with the 
known one-loop six-point next-to-MHV (NMHV) amplitude [33 133 UHl HI] • This agreement 
provides non-trivial evidence that there does indeed exist a Wilson-loop quantity dual to 
all scattering amplitudes. 

The aim of this paper is to combine various approaches in order to determine the six- 
point ratio function P, or equivalently the NMHV amplitude, analytically at two loops. This 
quantity was expressed in terms of dual conformal integrals, and computed numerically, in 
ref. [12]. We proceed in a manner similar to our recent examination of the three-loop six- 
point remainder function [22]. In particular, we make an ansatz for the symbols of the 
various pure functions involved. (See appendix [A] for a brief introduction to pure functions 
and their symbols.) In other words, we assume that the functions that appear fall within a 
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particular class of multi-dimensional iterated integrals, or generalized polylogarithms. We 
say functions rather than function because in general, beyond one loop, one can imagine that 
there are both (parity even) x (parity even) and (parity odd) x (parity odd) contributions, 
in a sense which we make specific in the next section. For convenience, we call these 
contributions 'even' and 'odd', respectively. At tree level and at one loop, the odd part 
vanishes. 

After constructing an ansatz, the next step is to impose consistency conditions. We first 
impose the spurious pole and collinear conditions. Then we impose that a certain double 
discontinuity is compatible with the OPE [23] |25| 155] . At this stage, we find that the 
symbols for the relevant functions contain nine unfixed parameters. We convert these sym- 
bols into explicit functions. In general, this step leads to 'beyond-the-symbol ambiguities'. 
These ambiguities are associated with functions whose symbols vanish identically, namely 
transcendental constants, such as the Riemann ( values (p, multiplied by pure functions 
of lower degree. However, in the present case, after re-imposing the spurious and collinear 
restrictions at the level of functions, there is only one additional ambiguity, associated with 
adding the product of (2 with the one-loop ratio function. This term obeys all constraints by 
itself and has vanishing symbol. We are thus left with a ten-dimensional space of functions. 
In particular, we find that the odd part is necessarily non-zero. Moreover, it is uniquely 
determined in terms of the even part. 

In order to fix the remaining free parameters, we turn to a representation of the even part 
of the two- loop six-point NMHV amphtude based on loop integrals [H]. We analyse this 
representation, appropriately rewritten with a mass regulator pT] , in the symmetric regime 
with all three cross-ratios equal to u. In this regime, the most cumbersome double-pentagon 
integrals can be traded for the MHV remainder function, plus simpler integrals. This 
observation allows us to perform an analytic expansion for small and large u. Comparing 
these expansions with the ansatz, we are able to match them, precisely fixing all remaining 
free parameters. The fact that the ansatz agrees with the expansion of the loop-integral 
calculation in this regime is a highly non-trivial cross check, since an entire function is 
matched by an ansatz with just a few free parameters. Further confirmation that our 
result is correct comes from comparing with a numerical evaluation [l2] at a particular 
asymmetric kinematical point. This latter check also confirms the expectation that V is 
defined independently of any infrared regularization scheme. See also ref. [25j for a recent 
discussion of different infrared regularizations and regularization-scheme independence. 

In contrast to the the two-loop six-point MHV amplitude |TB] , the two- loop six-point ra- 
tio function cannot quite be expressed in terms of classical polylogarithms. Two additional 
functions appear, one in the even part and one in the odd part. However, these functions 
have a very simple structure: we can write them as simple one-dimensional integrals over 
classical polylogarithmic functions of degree three. The even part of the ratio function 
can be written in terms of single-variable polylogarithmic functions whose arguments are 
rational in the three cross-ratios u, v, w, plus one of the new functions, which coincides with 
the finite double-pentagon integral fi*^^-* [57]. We use the differential equations obeyed by 
this integral [58] 159] to derive various parametric integral representations for it. The odd 
part consists entirely of the second new function, V, which also can be expressed as a single 
integral over classical polylogarithms of degree three. This function can also be identified 
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as the odd part of another finite double-pentagon integral, fi*^^^ [57], which we compute 
using the differential equations derived in ref. [58] . 

This paper is organized as follows. In section [2] we review non-MHV amplitudes, and 
the definition of the ratio function in planar = 4 super Yang-Mills theory. We discuss 
the physical constraints satisfied by P, namely in the collinear and spurious limits, and also 
those arising from the OPE expansion of super Wilson loops. We make an ansatz for the 
symbol of V at two loops in section [3l and then apply the constraints. In order to promote 
the symbol to a function, we introduce in section Hltwo new functions that are not expressible 
in terms of classical polylogarithms, but have simple parametric integral definitions. Next, 
in section [5l we parametrise the beyond-the-symbol ambiguities and apply the collinear and 
spurious constraints at the functional level, which leaves only ten unfixed parameters. In 
section E] we determine these parameters by performing an analytic two-loop evaluation of 
the integrals contributing to the even part of the NMHV amplitude in a special kinematical 
regime. The final result for the full two-loop NMHV ratio function is presented in section 
[71 We conclude in section [HI Several appendices contain background material and technical 
details. We provide the symbols for several of the quantities appearing in this article as 
auxiliary material. 



2 Non-MHV amplitudes and the ratio function 

To describe the scattering amplitudes of A/" = 4 super Yang-Mills theory, it is useful to 
introduce an on-shell superspace (see e.g. refs. [001 EB [3S1 [S2])- All the different on- 
shell states of the theory can be arranged into an on-shell superfield $ which depends on 
Grassmann variables 77"^ transforming in the fundamental representation of su{4), 

$ = 6-+ + t^^Va + ^V^'v^'Sab + Ir^^^'^^^'eABCDr'' + ^r^^r^^r/^r^^eABCDG-. (2.1) 

Here T a, Sab = ^^abcdS'~^^ , F"^, and G~ are the positive-hehcity gluon, gluino, 
scalar, anti-gluino, and negative-helicity gluon states, respectively. These on-shell states 
carry a definite null momentum, which can be written in terms of two commuting spinors, 
poa _ \ot^a^ Note that the spinors A and A are not uniquely defined, given p; they can 
be rescaled by A — j- cA, A — )■ c~^A. The transformation properties of the states and the 77 
variables are such that the full superfield has weight 1 under the following operator. 
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All the different (colour-ordered) scattering amplitudes of the theory are then combined 
into a single superamplitude A{^i, $2, • • • , $n), from which individual components can be 
extracted by expanding in the Grassmann variables rjf associated to the different parti- 
cles. The tree-level MHV superamplitude is the simplest cyclically invariant quantity with 
the correct scaling behaviour for each particle that manifests translation invariance and 
supersymmetry, 

A^'^ (2 3) 
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The arguments of the delta functions are the total momentum = A"Af and total 
chiral supercharge g""^ = ^tvf: respectively. The full MHV superamplitude is the 
tree-level one multiplied by an infrared- divergent factor, 

^MHV = ^Shv X M . (2.4) 



Moving beyond MHV amplitudes, we define the ratio function by factoring out the MHV 
superamplitude from the full superamplitude [3B], 

A = Amhy X V. (2.5) 

Here V has an expansion in terms of increasing Grassmann degree, corresponding to the 
type of amplitudes (MHV, NMHV, N^MHV, etc.), 

V = 1 + Pnmhv + Vn^muv + • • • + '^MHv • (2-6) 

The number of terms in the above expansion of V is {n — 3), where n is the number of 
external legs. The Grassmann degrees of the terms are 0, 4, 8, ... , (4n — 16). At six points, 
which is the case of interest for this paper, there are just three terms, corresponding to 
MHV, NMHV and N^MHV. The N^MHV amplitudes for n = 6 are equivalent to MHV 
amplitudes, which are simply related to the MHV amplitudes by parity. Thus the non- 
trivial content of the ratio function at six points is in the NMHV term. 

At tree level, V is given by a sum over dual superconformal 'i?- invariants' [36]. In 
particular, for six points we have 

"^NMHV ~ -Rl;35 + -Rl;36 + -Rl;46 ■ (2-7) 

The i?- invariants can be described using dual coordinates Xi, 6i defined by 

Pi - W - - Qi - \Vi - - ■ (2.8) 

Then we have [301 SS] 



(a. 


a-l){b,b 


-1) 6^{{r\ 




\dbr) - 


1- {r\XrbXba\ 


\6 


ar)) 






\-^ra-^ab\ 


\b- 


1) (^|^rb*^ba 






a 


\a-l) 



(2.9) 



The i?- invariants take an even simpler form in terms of momentum twistors [63l H5] . 
These variables are (super)twistors associated to the dual space with coordinates x, 9. They 
are defined by 

Z, = {Z,\x^), Zf="'" = (Ar,xfA.^), xf = ^r^A.„. (2.10) 

The momentum (super)twistors Zi transform linearly under dual (super) conformal sym- 
metry, so that {abed) = e^sTU Z]: is a dual conformal invariant. The i?-invariants 
can then be written in terms of the following structures: 

[abode] = jXajbcde) + cyclic) 

(abed) (bade) {cdea){deab){eabc) 
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which contain five terms in the sum over cychc permutations of a, b, c, d, e in the delta func- 
tion. The bracket notation serves to make clear the totally anti-symmetrised dependence on 
five momentum supertwistors. The quantity Rr-ab is a special case of this general invariant, 

Rr;ab = [r, o - 1, a, 6 - 1, b] . (2.12) 

At the six-point level it is clear that there are six different such invariants. We label them 
compactly by (t), using the momentum twistor t that is absent from the five arguments in 
the brackets: 

(1) = [23456], (2.13) 

and so on. 

In general i?- invariants obey many identities; see for example refs. [511 SI]- These 
identities can be organised as residue theorems in the Grassmannian interpretation [H]. At 
six points, the only identity we need is [36] 

(1) - (2) + (3) - (4) + (5) - (6) = 0. (2.14) 

Using eqs. fl212D . (|2J3|) and (ITO]) . we can rewrite the NMHV tree amphtude as 

^NMHV = [12345] + [12356] + [13456] = (6) + (4) + (2) = (1) + (3) + (5). (2.15) 



Beyond tree level, the i?-invariants in the ratio function are dressed by non-trivial func- 
tions of the dual conformal invariants [36]. In the six-point case, there are three independent 
invariants. We may parametrise the invariants by the cross-ratios, 

2 2 2 2 2 2 

u = ^f^, v = ^f^, w = ^f^. (2.16) 

•^U-^m •^Z&-^2b 

Often it will also be useful to use the variables yuiVv^yw defined by, 

U-Z+ V-Z+ w- z+ 

Vu = , Vv = , Vw = , (2.17) 

u — Z_ V — Z_ W — 



where 



^ 2 



A = {l-u-v-wy -Auvw. (2.18) 



In terms of momentum twistors, the cross-ratios are expressed as 

(6123)(3456) (1234)(4561) (2345)(5612) 

U = -, V = -, W = , (2.19) 

(6134)(2356) (1245)(3461) (2356)(4512) ' ^ ^ 

while the y variables simplify to 

_ (1345) (2456) (1236) _ (1235) (2346) (1456) _ (2345)(1356)(1246) 

^" ~ (1235) (3456) (1246)' ~ (1234) (2456) (1356) ' ~ (1345) (2346) (1256) ' 

(2.20) 
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In this form, it is clear that a cychc rotation by one unit Zi — \ Zi^i maps the y variables 
as follows, 

111 

Vu — > —, Vv — > — , Vw — > — , (2.21) 

Vv Vw Vu 

while the cross-ratios behave in the following way, 

u — > V, V — > w, w — > u . (2.22) 

The parity operation which swaps the sign of the square root of A {i.e. inverts the y 
variables) is equivalent to a rotation by three units in momentum twistor language. Indeed 
one can think of the cross-ratios as independent, parity- invariant combinations of the y 
variables. Specifically we have 

{1 - yuyv)il - VuVw) ' {I - yuy^){l - yuVw) ' 

and similar relations obtained by cyclic rotation. Because of the ambiguity associated with 
the sign of the square root of A in eq. 02.171) . the primary definition of the y variables is 
through the momentum twistors and eq. (12.201) . Further relations between these variables 
are provided in appendix |Fl 

At six points it can also be convenient to simplify the momentum-twistor four-brackets 
by introducing [16j antisymmetric two-brackets of CP variables Wi via 

(u) = iUjkimn{klmn) , (2.24) 

so that we have 

(12)(45) (24)(15) (26)(13)(45) 

(14)(25)' (14)(25)' ^" (46)(12)(35) ' ^ ^ 

plus six more relations obtained by cyclic permutations^ 

Having specified our notation for the invariants we need, we now parametrise the six- 
point NMHV ratio function in the following way, 

1 



Pi 



NMHV 



[(l) + (4)]V3 + [(2) + (5)]Vi + [(3) + (6)]r2 

+ [(1) - (4)]V^3 - [(2) - (5)]\4 + [(3) - (6)]^^] . (2.26) 

The Vi and Vi are functions of the conformal invariants and of the coupling, with the Vi 
even under parity while the Vi are odd (recall that parity is equivalent to a rotation by 
three units). The cyclic and reflection symmetries of the amplitude A (and hence the 
ratio function V) mean that the Vi and V are not all independent. Indeed, choosing 
V3 = V{u,v,w) and V3 = V{yu,yv,yw), we can write 

^ ^ [(1) + (4)]F(n, V, w) + [(2) + {5)]V{v, w, u) + [(3) + {6)]V{w, u, v) 



NMHV 



+ [(1) - mv{yu,y.,yn.) - [(2) - (5)]t^(y„ y^, + [(3) - (6)]V^(y^, 

(2.27) 



^In comparison with ref. [22|, the indexing of the Wi variables differs by one unit, and a square-root 
ambiguity in defining the y variables was resolved in the opposite way. 
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The functions V and V obey the symmetry properties, 

V{w,v,u) = V{u,v,w) , V{y^,y^,yu) = -V{yu,yv,yw) ■ (2.28) 

Note that we have written the parity-odd function V as a function of the y variables, while 
V, being parity even, can be written as a function of the cross-ratios. The functions V and 
V depend on the coupling. We expand them perturbatively as follows, 

oo oo 

V{a) = a'V'-'K ^(«) = Yl "'^^'^ • (2.29) 

1=0 1=0 

Here 

(2.30) 

where g is the Yang-Mills coupling constant for gauge group SU{N); the planar limit is 
— )■ oo with a held fixed. 

At tree level, we have V^^'' = 1 while V^(°^ vanishes. One can see that the expres- 
sion (I2.27P with V = 1 agrees with eq. fl2.15p . At one loop, V still vanishes, while V is a 
non-trivial function involving logarithms and dilogarithms, 

V^^^ = ^ -logulogtt; + log(uti?) logt; + Li2(l - n) + Li2(l - f) + Li2(l - w) - 2^2 , 

(2.31) 

V^^^ = 0. (2.32) 

The main results of this paper are analytical two-loop expressions for 1/(2) and V'-^\ both 
of which are non-vanishing. 

Let us discuss some general constraints that the functions V and V obey. 

Physical poles in amplitudes are associated with singular factors in the denominator 
involving sums of color-adjacent momenta, of the form [pi + pj+i + . . . + Pj-iY = ^Ij- Iii 
the /^-invariants, in the notation of eq. (12. lip , such poles appear as four brackets {abed) of 
the form 

{i - 1,1,3 - 1, j) = - (j - 1, j) 4 • (2-33) 
However, the i?-invariants also contain spurious poles, which arise from the four brackets 
{abed) that are not of this form. The full amplitude must not have such poles. Therefore 
the functions V and V must conspire to cancel the pole with a zero in the corresponding 
kinematical configuration. 

In the dual-coordinate notation (12.91) . the i?-invariants contain poles from denominator 
factors of the form {r\xraXab\b)- For special values of a,b,r, such factors can simplify into 
physical singularities, but for generic values they correspond to spurious poles. In the 
six-point case, for example, Ri-^q contains a factor of 

(l|xi4X46|5) = (l|xi4|4](45) (2.34) 

in the denominator. While the pole at (45) = is a physical (collinear) singularity, the 
pole at (l|xi4|4] = is spurious. In momentum-twistor notation, the spurious pole comes 
from any four-bracket in the denominator which is not of the form {i — — 1, j). 
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For example, in the six-point case the i?-invariants (1) and (3) both contain the spurious 
factor (2456) in the denominator. (In the dual-coordinate notation, this particular pole is 
proportional to (2|x25|5] rather than (l|a;i4|4].) In the tree-level amplitude fl2.15p there is 
a cancellation between the two terms, so we see that 

(1) ^ -(3) as (2456) ^ 0. (2.35) 

At loop level, using this relation, we find that the absence of the spurious pole implies the 
following condition on V and V, 

[V{u, V, w) - V{w, M, v) + V{yu, Vv, Vw) - V{yy,, yu, Vv)] (2456)=o = . (2.36) 

As the spurious bracket (2456) vanishes, we find the following limiting behaviour, 

,u(l — v) 1 (u — v)^ 1 — u 
w^l, yu^{l-w}- r^, Vv^T-. r, Vw ^ z • (2.37) 

This is easiest to see in the two-bracket notation of eq. f l2.24p . in which (2456) = cor- 
responds to (13) = and hence to Wt_= w^. The above condition reduces to the one of 
ref. |19] on the assumption that V = Oo The one- loop expression for V, eq. (I2.3ip . satisfies 
the above constraint with V^*^^) = 0, since logw — )■ in the limit. 

There is also a constraint from the collinear behaviour. There are two types of coUinear 
limits, a '/c-preserving' one where N^MHV superamplitudes are related to N'^MHV superam- 
plitudes with one fewer leg, and a '^-decreasing' one which relates N*^MHV superamplitudes 
to N^~^MHV superamplitudes with one fewer leg. These two operations are related to each 
other by parity and correspond to a supersymmetrisation of the two splitting functions 
found when analysing pure gluon amplitudes [511 E51 EZl ESI 112] • For the six-point NMHV 
case, we only need to examine one of the collinear limits; the other will follow automatically 
by parity. 

Under the collinear limit, the n-point amplitude should reduce to the {n — l)-point 
one multiplied by certain splitting functions. The splitting functions are automatically 
taken care of by the MHV prefactor in eq. (12. 5p . The n-point ratio function V should 
then be smoothly related to the {n — l)-point one. Consequently, in the collinear limit 
the loop corrections to the six-point ratio function should vanish, because the five-point 
ratio function (containing only MHV and MHV components) is exactly equal to its tree- 
level value. The i?-invariants behave smoothly in the limit, either vanishing or reducing 
to lower-point invariants. In the case at hand we can consider the limit Zq ^ Zi, which 
also corresponds to Wq — )■ Wi, or x^^ — )■ 0, or w — )■ with v ^ 1 — u. In this limit, all 
i?- invariants vanish except for (6) and (1), which become equal. Beyond tree level, the 
sum of their coefficients must therefore vanish in the collinear regime. This implies the 
constraint, 

[V{u, V, w) + V{w, u, v) + V{yu, yv, yw) - V{yy,, yv)]w^Q, v^i~u = . (2.38) 

In fact the parity-odd function V drops out of this constraint. The reason is that the 
collinear regime can be approached from the surface A{u, v,w) = (see eq. f l2.18p ). and all 
parity-odd functions should vanish on this surface. 



^This is true after correcting a typo in eq. (3.63) of that reference. 
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The final constraint we will need comes from the predicted OPE behaviour of the ratio 
function [55]. The general philosophy that an operator product expansion governs the 
form of the amplitudes comes from the relation of amplitudes to light-like Wilson loops. 
Light-like Wilson loops can be expanded around a coUinear limit and the fluctuations can 
be described by operator insertions inside the Wilson loop [221 ES]- By extending 
this philosophy [55] to supersymmetrised Wilson loops [511 [50] (or equivalently correlation 
functions [52]) one can avoid questions about giving a precise Lagrangian description of 
the object under study. In this sense the OPE can be used to justify the existence of a 
supersymmetrised object dual to non-MHV amplitudes. 

The analysis of ref. [55J allows one to choose various components of the ratio function. 
Let us consider the component proportional to X2X3X5X6- The only term in eq. f l2.27p that 
contributes to this component is the first one, 

^Si= (2.39) 

In order to examine the OPE, we follow ref. [55] and choose coordinates (r, a, 0) by fixing 
a conformal frame where 



1 1 \fmyw 



(2.40) 



(2356) 4(cosh 0" cosh r + cos 0) 2(1 — t>) 
and the three cross-ratios are given by 

e'^ sinh r tanh r 1 e"*^ sinh r tanh r 

XL = V = VJ = 

2 (cosh o cosh r -|- cos 0) ' cosh^ r ' 2 (cosh a cosh r + cos 



Extrapolating the results of ref. [55j to two loops, the OPE predicts the leading (double) 
discontinuity of the (2356) component of the ratio function to be. 



m=— OO 



where 

^tT^ (^) = sech^r 2F1 [i {E - ip) , |(E + ip) ; E; sech V] , (2.43) 

4''''^^) = \{-^TB [i(|m| + 1 + ^p), |(|m| + 1 - , (2.44) 

7i+|m|(p) = V^(i(l + |m| +ip)) +V^(^i(l + |m| -zp)) -2V^(1) . (2.45) 

Here 2-^1 is the hypergeometric function, i?(a,/3) = r(a)r(/3)/r(a + /3) is the Euler beta 
function, and t\} is the logarithmic derivative of the F function. 



3 Ansatz for the symbol of the two-loop ratio function 

In order to make a plausible ansatz for the ratio function at two loops we assume that 
the functions V'^'^^iu^v^w) and V'^'^\u^v^vS) are pure functions of m, v and w, i.e. iterated 
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integrals or multi-dimensional polylogarithms of degree four. Moreover, we make an ansatz 
for the symbols of V^'-^^ and V^'^\ requiring that their entries are drawn from the following 
set of nine elements, 

{u,v,w, l-u,l-v,l- w,yu,yv,yw} ■ (3.1) 

We summarise some background material on pure functions and symbols in appendix El 
We recall that the y variables invert under parity. The parity-even function V should have 
a symbol which contains only terms with an even number of y entries. Likewise the symbol 
of the parity-odd function V should contain only terms with an odd number of y entries. 
The ansatz for the symbol entries is the same as the one used recently for the three-loop 
remainder function [22J (after omitting restrictions on the final entry of the symbol). It 
is consistent with every known function appearing in the six-point amplitudes of planar 
A/" = 4 super Yang-Mills theory, in particular the simple analytic form of the two-loop 
remainder function found in ref. [16]. It is also consistent with the results for explicitly 
known loop integrals appearing in such amplitudes, see refs. [58l [59l [67] . In the ensuing 
analysis we will find many strong consistency checks on our ansatz. 

Let us pause to note that our assumption that the relevant functions are pure functions of 
a particular degree equal to twice the loop order is by no means an innocent one. Although 
it is true that such general polylogarithmic functions generically show up in amplitudes in 
four- dimensional quantum field theories, it is certainly not true in general that they always 
appear with a uniform degree dependent on the loop order. In QCD, for example, the 
degrees appearing range from twice the loop order to zero, and the transcendental functions 
typically appear with non-trivial algebraic prefactors. In fact the observed behaviour of 
having maximal degree only is limited to A/" = 4 super Yang-Mills theory, and the most 
evidence is for the planar sector. This behavior is the generalization, to non-trivial functions 
of the kinematics, of the maximal degree of transcendentality for harmonic sums that has 
been observed in the anomalous dimensions of gauge- invariant local operators [68] . 

The symbols we construct from the set of letters fl3.ip should obey certain restrictions. 
They should be integrable; that is, they should actually be symbols of functions. The initial 
entries of the symbol should be drawn only from the set {u, v, w}, because the leading entry 
determines the locations of branch points of the function in question, and branch integrable 
symbols of degree 4 for V, and 2 for V, obeying the initial entry condition as well as the 
symmetry conditions fl2.28p . The spurious pole conditions fl2.36p provide 14 constraints and 
the coUinear conditions (12.380 provide 14 more, leaving 15 free parameters at this stage. In 
order to impose the constraints from the leading discontinuity predicted by the OPE, we 
use the fact that the sum (12.420 is annihilated by the following differential operator [55] , 

V = d^ + 2coth{2T)dr + sedi'^Tdl + dl + l. (3.2) 

In the u,v,w variables this differential operator is given by 

V = |(I?+ + I?_) + 1, (3.3) 
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where 



1-v 



— z±udu — (1 — v)vdv — z±wdw 

+ (1 — u)vuduudu + (1 — v^vd^vdv + (1 — w)vwdwwdw 

+ (—1 + u — V + w){{l — v)uduvdv — vuduwdw + (1 — v)vdvwdw) 



(3.4) 



Imposing that the double discontinuity f l2.39p is annihilated by the operator V gives 
5 further conditions, leaving 10 free parameters in the symbol. One of these parameters, 
denoted by ax below, is just the overall normalisation of the symbol of the double discon- 
tinuity, which is non-zero and convention-dependent. In the following we provide functions 
with physical branch cuts which represent the symbol. We find that the solution has the 
form, 

9 

S{V) = ax5(Vx) + ^a,5(/,), S{V) = ax S{Vx) + asS{f) , (3.5) 

i=l 

where ax and ai through ag are the constant free parameters, and the quantities V, fi, 
Vx and / will be defined belowjf] 

The double discontinuities of the functions appearing in eq. f l3.5p obey 

S{AyAvVx) = 2 ax u(8)(l — n) + u^u + w®{l — w) + w ®w 

+ 2(^u 1^ w + w 1^ u — uw ® (1 — f ) — (1 — f ) uw 

+ (l-t;)®(l-t;))] , (3.6) 

S{A,AJ\) = 0, (3.7) 
S{A,KVx) = SiA.Aj) = 0. (3.8) 



Consistency with the spurious pole condition fl2.36p forces the odd part V to be non-zero, 
given that ax is non-zero. The odd part contains no ambiguity at the level of the symbol 
(or beyond it), once we fix the even part, particularly the two parameters ax and as. 

The symbol of the double discontinuity of V and V is entirely controlled by Vx , through 
eq. (13. 6p . We can find a function compatible with this symbol, and compare it to eq. (I2.42p 
to fix the (2 terms. We find that 



'-^v'-^v' NMHV ^ 



1 



(2356) 



log^ u + log^ w + 4: log u log w + 2 log^(l — v) 



- 4 \og{uw) log(l -v)-2 (^Li2(l -u) + Li2(l - w) - 2 C2 



(3.9) 



In order to present the symbols appearing in eq. (13. 5p explicitly and compactly, it is very 
useful to employ harmonic polylogarithms [70l [7T| [72] . This presentation simultaneously 

^In section [6] we will fix the ten parameters using an analytical computation for particular kinematics. 
We have compared the symbols p.Sp for V and V with all parameters fixed to an independent computation 
of these symbols from a formulation of the super Wilson loop ^ ; the results agree precisely. 
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accomplishes the following step, of turning the symbols into functions, up to certain beyond- 
the-symbol ambiguities. The functions we will present are of degree at most four, and almost 
all of them can be represented in terms of classical Li„ functions. Thus the use of harmonic 
polylogarithms may seem unnecessarily complicated. However, it is a very useful way to 
represent, at any degree, a symbol only involving the letters {u, v,w,l — u,l — v,l — w}, 
whereas Li„ functions are often insufficient beyond degree four. 

Harmonic polylogarithms are single-variable functions defined by iterated integration. 
It is very simple to write down their symbols. We use harmonic polylogarithms with labels 
(weight-vector entries) "0" and "1" only. The symbol of a harmonic polylogarithm of 
argument x is obtained by reversing the list of labels and replacing all "0" entries by x and 
all "1" entries hy 1 — x. Finally, one multiplies by (—1)" where n is the number of "1" 
entries. For example, the symbol of i?o,o,o,i(3^) — Li4(a^) is — {1 — x) <S> x <S> x <S> x, and the 
symbol of i?o,i,o,i(3^) is {1 — x) <S> x <^ {1 — x) <^ x . 

We also use the common convention of shortening the label list by deleting each "0" 
entry, while increasing by one the value of the first non-zero entry to its right, so that, for 
example, i7o,o,o,i(2;) = H^i^x) = Li4(x) and ifo,o,i,i(^) — -f^3,i(^)- Apart from the logarithm 
function, we take all arguments of the harmonic polylogarithms to be (1 — -u), (1 — v) 
or (1 — w). This representation guarantees that the functions we are using to represent 
the symbol do not have any branch cut originating from an unphysical point. We then 
compactify the notation further by writing i?2,2(l — x)— H2 2, and so on. Finally, we recall 
that the symbol of a product of two functions is given by the shuffle product of the two 
symbols. 

With this notation we can immediately write down a function which has the symbol, 
S{Vx), of the part of V with non-zero double discontinuity, i.e. the part fixed by the OPE, 

Vx = {^Hl, + log u {H; + 2Hl^ - + QH^^ + fi/^" log w) + log^ u {H^ - 3H^) 

+ log V [H^ - SH^, - I log « {H^ + HI) - I \oi u log iv] + log^ v {-H^ + i log^ u) 
+ {u'(r^w)^+ 4if^ 1 + 2\ogu log^ vlogw-^ log^ u log^ w . (3.10) 

We can similarly write down functions with the correct symbols for the first seven ambigu- 
ities in the even part (the double v discontinuity of each function vanishes): 

/TJU TJW 

/2 = [ - \ogu{H^ + H^^ + Hl^Xogw) - log^uiH^ + i log vlog^/;) + {u ^ w)] , 
/a = [ — H2 log ulogv + {u -ir^ w)] + H2 log u log w , 

f4—[ — H2 log U log W — log^ U {2H2 + log V log w) + {u -^r^ w)] — H2 log U log W , 
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/5 = [H^H^ + Hl^ + \ogu{2Hl^ - - \ogw) + \ogv {2Hl^ + \ogu {H^ + H^,)) 

+ {u^ w)] + , 

h = [- 2H^H^, - - + \ogu {-2H^ - 2HI, + 2H^, + H^logw) 

- \ogv {2HI, + \ogu (H^ + HI)) + {u^ w)] - 2HI, - AHl, - 2HI \ogv , 



h = [- - + log u {HI - 2HI, + H^^ + Hl^ log w) + log^ u {-^H^ + H^ 

+ log V {H^ + \ log2 u\ogw)^[u^w)\- mi - mi 1 1 - 2HI ^\ogv - \Hl log 



For the even part there remain two more ambiguities whose symbols cannot be expressed 
in terms of those of the single- variable harmonic polylogarithms with the arguments we have 
been using, 

h=\- HlHl - 2HI + Hl^ - AHl, + QHl,^, + log u {-HJ^ - HI + H^ + 2HI,) 
+ Hl^ log^ M + log {H^ - H^ log u) + {u ^ w)] 
-Hl^ - 2Hl^ + H:^ log M log w + I log^ u log^ w-Hl\ogv- 2^](2) (^^ v) , 

h = T^'i\u,v,w) . 

Here TZ^ stands for the two- loop remainder function, whose symbol is known |T6] . The 



appearance of the two-loop remainder function as an ambiguity should not be surprising. 
It is a function with physical branch cuts, which vanishes in the coUinear limit. Also, it 
is totally cyclic and hence automatically satisfies the spurious pole condition on its own. 
Furthermore, it has vanishing double discontinuities and hence drops out from the leading- 
discontinuity OPE criterion (12.421) . It is known [16] that T^g^'' can in fact be expressed in 
terms of single- variable classical polylogarithms. However, to do so one must use arguments 
involving square roots of polynomials of the cross-ratios. 

The other quantity not given in terms of the harmonic polylogarithms, which enters fs, is 
the integral In fact its symbol can also be recognised from other considerations |58l[59] . 
as we will discuss in the next section. The symbol of fi'-^^ is. 



2 

where 

1111 II I ^ — 1)1 

(3.12 



5(fi(2) ^u, V, w)) = -I S{q^) ® + S{qr) ® r + S{^q) ® VuVv , (3.11) 

uv u{l — v) 

~ {l-u){l-v)' ^ ~ v{l-u)' 

Here $6 is the one- loop six-dimensional hexagon function [591 EZ], whose symbol is given 
explicitly in terms of the letters of our ansatz |59J, 

S{^g) = -S{Q'-^\u,v,w)) ®?/^ + cyclic, (3.13) 

where Q,^^^ is a finite, four-dimensional one-loop hexagon integral [STj 158]. 

n^-^\u,v,w) = log M log + Li2(l -u) + Li2(l - v) + Li2(l -w)- 2(2 . (3.14) 
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The other degree 3 symbols above can be represented by harmonic polylogarithms as follows, 
= [-H^ - - HI log M - i log' u\ogv + log w + {u ^ v)] 
+ 2H^-^ + hgw + logM logt; log w , 
Qr = [-H^ + Hl^ + H^logu + H^logu + llog^ulogv-iu^v)]. (3.15) 

For the symbol of the parity-odd function V, we find that the part fixed by the OPE 
(acting in conjunction with the spurious-pole constraint fl2.36p ). SiVx), coincides with the 
symbol of 

Vx = l-elogf-). (3.16) 
For the odd part of the ambiguity associated with «§, we have 

I — u 

S{f) = S{fu)^yu + S{U)^y, + S{U)^y^-S{%)^- , (3.17) 

1 — w 

where the functions fu, fv, fw are given by, 

A = [2H^ - \ogv-{ui^w)]- 2HI, - H^, log v , 

A, = [2H^ - 2Hl^ - log u - HI log u + Hl^ log u - [u ^ w)] , 

U = [2H^ - H^ log v-{u^w)] + 2Hl^ + H:^ log v . (3.18) 

We emphasise again that the formulas presented in this section are meant to represent 
the symbols of the functions involved. For some of the relevant symbols [SiVx) and S{fi) 
through S{ff)) we were able to trivially write down actual functions which represent those 
symbols in terms of single- variable harmonic polylogarithms with arguments 1 — x, where 
X is one of the cross-ratios. For two others {S{fg) and SiVx)) we recognised them as 
involving symbols of functions we already know, namely the two-loop remainder function 
and the one-loop six- dimensional hexagon integral. In order to write down actual functions 
for V and V there are two issues to address. Firstly, we must give functions which represent 
the symbols S{n^'^\w,u,v)) and S{f). Secondly, we must include all possible terms which 
have vanishing symbol and which are therefore insensitive to the analysis we have presented 
so far. We address these two issues in the next two sections. 

4 Digression on integral representations for and / 

Here we will present integral formulas to define the functions i7^'^ and / whose symbols 
are given in the previous section. We also present a new representation of the two-loop 
remainder function, based on the integral This section is more technical, and could 

therefore be skipped on a first reading. 

4.1 Integral representations for il^^^ 

We start with the finite double-pentagon integral ^l^'^^u, v, w) [57J. Let us take a derivative 
with respect to w. The only contributing term from the symbol fl3.1ip is the last one, so 
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we see that the symbol iS(n*^^^) is consistent with the differential equation, 



> (m, v,w) = — —d^ logiyuVv) 



(4.1) 



We recognise here the differential equation |59j relating the two-loop, finite double-pentagon 
integral to the massless, one-loop, six- dimensional hexagon function $6. 



The relation (14. ip can be used to write an integral formula for 

dt 



^/A{u,v,t) 



(4.2) 



The relevant boundary condition is il^'^\u,v,0) = '^^'^\u,v), where \E'*-^^ is the two-loop 
pentaladder function found in ref. [58]. The boundary behaviour at w = was tested 
numerically from the Mellin-Barnes representation for fi*^^) [58]. The symbol (13. lip reduces 
to the symbol of at w = 0. It is the unique symbol within our ansatz, built from the 
letters in eq. (13. ip . that obeys eq. (14. ip and the w = boundary condition. The integral in 
eq. (14.20 is well-defined and real in the Euclidean region, i.e. the positive octant in which 
u,v,w are all positive, because the integrand $6 = ^e/v^ is well-defined and real there, 
and $6 is well-behaved even where A vanishes [59]. 

As discussed in ref. [SSj , the first-order differential equation (14. ip can be obtained from 
the second-order equation of ref. [58] for the double-pentagon integral, which can be written 

as 



-u 



(1 — u)du — v{l — v)d^ + {1 — u — v){l — w)dyj Vl^'^\u., w) {u,v,w) . (4.3) 



Because the second-order operator naturally factorises into two first-order operators, we 
can integrate up to fi*-^-* in two steps. This procedure will yield another one-dimensional 
integral relation for fi^^^. We define 

Q^{u,v,w)= —u{l~u)du — v{l — v)dti + {l—u — v){l — w)dwQ^'^\u,v,w), (4.4) 

so that 



The above formula can be used to define the function Q^, 

Li3(l -w) + Li3 ( 1 - - 

w 



(4.5) 



(4.6) 



Li2(l - w) + Li2(l - u) + Li2(l - v) + \ogu logv -2(2 

1 

u 

112(1 -u) - 112(1 - V 



+ log w 



--log3tf;-2Li3(l-n)-Li3( 1 



2Li,(l -1;) - Li, 1 



+ log - 

logM logf log(Mf) . 

2 



+ - log M + - log V 
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This function obeys eq. (14 .Sp and has a symbol coinciding with that of from eq. (I3.15p . 
It also obeys Q(f,{l, 1, 1) = 0. In principle, eq. (14. 5 p allows one to add beyond-the-symbol 
terms to that are proportional to (2\og{uv), and to (3. We verified numerically that 
these terms are absent. The function is manifestly real in the positive octant. 

Given the function Qtf,, we can integrate eq. (14. 4p to obtain fi*^^). We first note that the 
relevant operator becomes very simple in the {yu,yv,yw) variables, 

- u{i - u)du - v{i - v)d, + (1 - « - ^;)(l - w)d^ = - y^)i^ - y^y^y^) g (4.7) 

Inserting this relation into eq. (14. 4p . we find an alternative integral formula for ^l^"^^ in terms 
of the y variables, 

n('\u,v,w) = -6C4+ r -^^^^^Myu.y..t) . (4.8) 

1 - 1 1 - yuyvt 



yuvv 



Here we use the notation Q^(?/„, yjj = Q<f>{u{yu, yv, yw), v{yu, yv, yw), w{yu, yv, yw))- Note 
that this integral is well-defined at the lower limit of integration, for the following reason: 
Whenever the product of the y variables is unity, Vuyvyw = 1, we see from eq. (I2.23P that 
the cross-ratios collapse to the point (M,f,w) = (1,1,1), and at that point vanishes, 
(5,^(1, 1, 1) = 0. Equation (14. 8 p can be applied straightforwardly in the y variables for 
yw > 1 and < 1, and also for y^ < 1 and yuyv > 1- (In other regions, the vicinity of 
t = 1 makes a direct integration problematic.) 

It can be more convenient to map the integral (14. 8 p back to the (m, v, w) space. This 
mapping avoids problems related to the variables {yu,yv,yw) becoming complex when A is 
negative. To do this mapping, we first define 

u{l-v) yuil-yvf 



v{l-u) y^il-yuY' 
u{l -u)v{l-v) ^ yujl -yufyvil - yvf 
{i-wf {i-yuyvY 



(4.9) 



(4.10) 



^ ^ 1 - w ^ (1 - yuyv) (1 - yuyvyw) ,^ 
uv yu{'^ - yu) yv{'^ - yv) {'^ - yw) ' 

Notice that r{yu,yv,yw) and s{yu,yv,yw) are actually independent of yw Therefore the 
curve of integration in the integral (14. 8 p from (1,1,1) to {u,v,w), which has constant yu 
and yy, should have a constant value of r and s along it, while t varies. Also, 

dOogt) ^ 1 - Vuyv ^2) 

dVw (1 - y«;)(l - VuVvyw) ' 
so that the measure in eq. (14.81) is just dlogt. 

Let {ut,Vt,Wt) be the values of {u,v,w) along the curve from (1,1,1) to {u,v,w). We 
solve the two constraints, that r and s are constant along the curve, i.e. 

ut(l-vt) u(l-v) 



Vt{l - ut) v{l - u) ' 

ut{l - Ut)vt{l - vt) u{l-u)v{l 



Wt)'^ (1 — w)" 



(4.13) 
(4.14) 
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for Vt and Wt in terms of Ut, obtaining, 



{l-u)vut 

Vt = —J. ^ — , 4.15 

u\ \ — v) ^ \ v — U) Ut 

{l-w)ut{l - Ut) 
Wt = I T- ^ . 4.16 

U [1 — V) + [V — U) Ut 



Inserting these expressions into cilogt = (ilog[(l — Wt)/ut/vt], we have 

d{\ogt) _ 1 
dut utiut - 1) ' 

which enables us to use Ut as the final integration parameter, 



(4.17) 



n^''\u,v,w) = -QU+ r —^^^Q^iut,vt,wt). (4.18) 

Jl Ut[Ut - 1) 

Using this formula, with Q^j, from eq. (14. 6p . for which the polylogarithms are all rational 
functions of the cross ratios, it is easy to rapidly get high-accuracy values for VP''^. For 
example, we find 

l](2)(ff,f ,^) = -5.273317108708980008, (4.19) 
fi(2)(f,^,ff) = -6.221018431345742955, (4.20) 
l](2)(^,ff,f ) = -9.962051212650647413, (4.21) 

in general agreement with the numbers obtained at these points using a Mellin-Barnes 
representation for the loop integral. 



4.2 A new representation of the two-loop remainder function 

Now that we have obtained representations of the function we note that the two- loop 
remainder function can be written in terms of this function, together with functions with 
purely rational (y-independent) symbols. Specifically, we have 



nf\u,v,w) = ^ ^ 



Vt^'^\u,v,w) + Vt^'^\v,w,u) + Vt^'^\w,u,v) 



The piece with a rational symbol is defined as 



(2) 



'6. rat 



^(Li2(l-l/M)+Li2(l-l/t;)+Li2(l-l/w))'+r(M)+r(t;)+r(«;)-C4 



with 



1 



(4.22) 



(4.23) 



riu) = — Li4(M) — 114(1 — u) + 114(1 — l/u) — logM 113(1 — l/u) — - log^w log(l — u) 

6 

+ ^ (Li2(l - V«)) ' + ^ log^ M + C2 (Li2 (1 -u)+ log2 + Cs log u . (4.24) 
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The function T^g l^^^ is real when all three cross-ratios are positive. Almost all of the terms 
in eqs. fl4.23p and f l4.24p make this manifest term-by-term, because they contain only loga- 
rithms of cross-ratios, or Li„(x) for some argument x which is less than one. The one slight 
exception is the combination 

Lu{u) + - log^ u log(l - u) . (4.25) 
6 

It is easy to see that eq. ( I4.25P is real as well, but in this case the branch cut starting at 
M = 1 in each term cancels in the sum. 

In one sense, the representation fl4.22p is a step backward from ref. [IS], because the 
function Q^'^\u,v,w) cannot be expressed in terms of classical polylogarithms, whereas 
TZ^Q ^ can be. (The absence of a classical polylogarithmic representation for Q'^'^^ can be 
seen from its symbol, using the test described in ref. [16].) However, the appearance of 
the sum over cyclic permutations of the finite two-loop double-pentagon integral is natural, 
and the coefficient of j matches the one in the expression for the two-loop MHV amplitude 

in ref. [57]. The relation fl4.22p between TZ^^ and fi*-^^ will be useful for us in the ensuing 
NMHV analysis. 



4.3 An integral representation for / 



We can obtain in a similar way an integral formula for the parity-odd function /. Note 
that we already have a formula, eq. (13.160 . for the function Vx = ^elog^u/w). It is useful 
to observe that the combination Vx + f has a symbol which can be arranged so that the 
final entries are drawn from the list, 



u{l — w) 



w{l — u) 

In terms of the y variables, the last final entry in the list above is independent of 

u{l - w) _ yu(l - ywf 
w{l-u) yu,{l-yuY' 



(4.26) 



(4.27) 



This fact allows us to obtain the symbol of the logarithmic derivative with respect to y^, 
which is independent of the y variables. 



S{Z) 



2{u 



u 



l-u 



:i-u)- 



u ® 



w 



l-u 



v®{l — v)+w®u 



® u 



+ {u®v + v®u) ® {1 — v) 



u u (^w 



(4.28) 



The combination Vx + f can then be written as an integral of a function with this symbol, 

Vx + f= I jZiy^,t,y^). (4.29) 
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Here Z is to be considered as a function of the variables y^) for the integration, but 

it is most simply expressed in terms of the variables 



Lis 1 - - 

u 



Lia 



1-- +log - 
w / \w 



u 



Li2(l-t;)-2Cs 



(4.30) 



in which form it is manifestly real in the positive octant. 



Using the same trick we used for we can rewrite the integral fl4.29p directly in the 
{u,v,w) space. The only difference is that the roles of {v,yy) and {w,yw) are swapped, and 
there is an extra factor multiplying the pure function, corresponding to 



(1 - yy){l - yuyyyw) 



(4.31) 



Thus we get. 



Vx + f 



dut 



ll Ut{Ut - 1 

-^/X{v~u~w) 



Z{ut,vt,wt) 



Vt 

dut Z(ut, Vt, Wt) 



where 



Vt 
Wt 



Vt [u {1- w) + {w - u) Ut\ 



[I - v)ut (1 - Ut) 



(4.32) 
(4.33) 



u{\ — w) + {w — u) Ut 

(1 — u)wut 
u{l — w) + {w — u) Ut 



(4.34) 
(4.35) 



The second form of the integral, eq. fl4.33p . makes clear that in the positive octant, Vx + f 
is real for A > 0, and pure imaginary for A < 0. The overall sign of eq. (I4.33P corresponds 
to the branch of a/A defined in term of the y variables in eq. (IF.SP : it ensures that the 
logarithmic derivative with respect to y^ reproduces Z. 



5 Ansatz and constraints at function level 

Now that we have obtained explicit functions representing the symbols in section [3l we 
proceed to enumerate the additional possible contributions, all of which have vanishing 
symbol. The ratio function is real-valued in the Euclidean region in which all three cross- 
ratios are positive. Each of the above functions entering V also has this property. Therefore 
any additional functions that we add to our ansatz must also obey the property. In addition 
to the parameters {ax, ai, • • • , ctg}, we have the following real-valued parity-even beyond- 
the-symbol ambiguities: 
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• At the (2 level, 

g^"^^ = (2 [ci (log^ u + log^ w) + C2 log^ f + C3 log('uw) log v + C4 log u log w 

+ c^H^ + H^) + CeH^2] , (5.1) 

• At the Cs level, 

g^^^ = Cs [c7 log(Mw;) + cg log t;] , (5.2) 

• At the ^4 level, 

^(^)=C4C9. (5.3) 

If our ansatz is correct, then we expect that the parity-even function V should be given 

by 

9 

1=1 

for some rational values of the and Cj. There are no parity-odd beyond-the-symbol 
ambiguities that possess only physical branch cuts. (This fact follows from the absence of 
an integrable parity-odd degree-two symbol whose first slot is constrained to be u, v or w.) 

Next we would like to apply the constraint fl2.38p from the coUinear limit, namely 
V{u,v,w) + V{w,u,v) — 7- as 1(7 — )■ 0, w — )■ 1 — u, but now at the level of functions, 
not just symbols. One way to do this is to first complete the functions Vx and fi into 
new functions Fx = Vx + Vx, Fi = fi + fi, each of which gives a vanishing contribution to 
V{u, V, w)+V{w, u, v) in the coUinear limit. Although the symbols of the functions /i, . . . , /g 
were already constrained to give a vanishing contribution in this limit, that does not mean 
that they vanish as functions. Instead we will correct Vx and the fi by appropriate beyond- 
the-symbol terms, Vx and fi, which are constructed from the expressions (15. ip . (15. 2p and 
(15. 3 p for suitable values of the constants Cj. The function /g requires no such correction, 
because it is the two-loop MHV remainder function, which vanishes in all coUinear limits. 

To perform this correction, we need to know the coUinear limits of the functions Vx 
and fi. For all but /§, these limits are straightforward to compute. The hmit of /§ is 
more complicated to obtain due to the presence of We compute this limit directly in 
appendix IHl However, we may also observe that in the coUinear constraint equation (I2.38p . 
only the combination ^2^^^ {w,u,l — u)+ fi^^^ (1 — u,w,u) is needed for smaU w. (To see this, 
we use the symmetry of fi*-^-* under exchange of its first two arguments.) This combination 
appears on the right-hand side of eq. (I4.22p . evaluated in the limit w 0, v ^ 1 —u, along 
with fi^^) {u, l — u,w) and the simpler function TZf l^^{u, l — u,w). Now the left-hand side of 
this equation, the two-loop remainder function, vanishes in the limit. Also, Q^'^\u, 1—u, 0) = 
\1/^^^(m, 1 — u) = 0, where '^^'^\u,v) is the two-loop pentaladder function [58]. Hence the 
limit of the pair of functions appearing in fsiu, v, w) + fs{w, u, v) reduces to evaluating 

(2) _______ 

^6 rat (^5 1 — M,0), usiug eq. (14.231) . The explicit formulas for all the required beyond-the- 
symbol functions Vx and fi are given in appendix [Bl 
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The collinear constraint at the level of functions fixes 7 out of the 9 beyond-the-symbol 
terms, leaving only the following combinations which have vanishing collinear contributions: 



~9i = C2[C2 + m-H^-Hn , (5.5) 
92 = C2 [-C2 + ^H^ + \og{uw) log f - log u log w] . (5.6) 

The function V should therefore be given by 

9 

V{u, V, w) = axFx + ^ aiFi + cigi + £2^2 , (5.7) 

i=l 

where ci and £2 are arbitrary constants. 

We can analyse the constraints coming from the spurious pole condition (I2.36P in a 
similar way. The end result of this analysis is that one more beyond-the-symbol ambiguity 
is fixed, leaving just one such function free. In fact, this is the maximum number of beyond- 
the-symbol terms we can fix with this analysis, because one can always add C,2 multiplied by 
the one- loop ratio function, V'^^\ given in eq. (12.3 ip . This product is a linear combination 
of gi and ^2, namely ^2 — Qi-, and it automatically satisfies all constraints by itself. Thus 
the only remaining beyond-the-symbol ambiguity is C2 V^^\ 

In the next section, we will calculate analytically the loop integrals contributing to 
the two-loop NMHV amplitude for special kinematics. We will use this information to 
determine the remaining unfixed parameters, ax-, through ag, and (one of) Ci and 62. 

6 Analytic calculation using loop integrals 

In this section, we will fix the remaining undetermined parameters in our ansatz by com- 
puting the ratio function analytically in a certain kinematical regime. 

We find it convenient to perform our calculation using a mass regulator [2T]. As was 
reviewed in section [2l the ratio function is infrared finite. Moreover, it should be indepen- 
dent of the regularization scheme used to compute it. We first verify this statement at one 
loop by re-evaluating the MHV and NMHV six-point amplitudes in the mass regulariza- 
tion. At two loops, we find agreement with previous numerical results |12] obtained using 
dimensional regularization. 

6.1 Review of six-point MHV amplitudes 

Recall from section [2] that supersymmetry allows to write any MHV amplitudes to all 
loop orders as a product of the tree-level amplitude, multiphed by a helicity-independent 
function. We have 

^MHv(a) = ^Siiv X M{a), (6.1) 
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where M{a) = 1 + aM^^^ + a^M^'^'^ + . . ., and a is defined in eq. fl2.30p . The known structure 
of infrared divergences takes a particularly simple form if we consider log M, namely [73] 



logM(a,4) 



E 



i=l 



2 -^1,1+2 



Qoia) 



log 



X. 



171^ 



16 m? 2 

Here 7(a) is the cusp anomalous dimension [73]. It is given by 

7(a) = 4a - 4(20^ + O(a^) , 

and we have 

Qo{a) = -C3a' + 0(a'), 
/(a) = ^a2 + C(a3). 



+ /(«) 



(6.2) 

(6.3) 
(6.4) 



Moreover, the finite part F satisfies a dual conformal Ward identity [U |5], whose most 
general solution is 



F(a,4) 



^7(a)FW(4) + 7^6(M, V, w; a) + C{a) + ^(m^ 



with 



C{a) = -^a^ + 0{a^). 



(6.5) 



(6.6) 



The first term on the right-hand side of eq. (16.51) comes from the BDS ansatz [12], and 
provides a particular solution to the Ward identity. It is given by the one-loop contribution 
to F, multiplied by one quarter of the (coupling-dependent) cusp anomalous dimension 7(a). 
Hence its kinematical dependence is determined by the one- loop result. The second term 
on the right-hand side of eq. (16. 5p . the remainder function TZQ{u,v,w;a) [10l[Tl], depends 
on three conformal cross-ratios, u, v and w. There is no remainder function for four and 
five points, because non- vanishing conformal cross-ratios only appear starting at six points. 
The specific choice of the kinematic-independent terms /(a) and C*(a), determined by the 
four-point and five-point cases, was made in such a way [12] that 71q{u, v, w; a) vanishes in 
the collinear limit. 



6.2 Six-point NMHV amplitudes and ratio function 

The NMHV amplitude can be written as 

^NMHv(a) = i^JSHvf[(2) + (5)]W^i(a) - [(2) - (5)] lyi(a) + cyclic 



(6.7) 



where Wi{a) = 1 + aW^^^ + a^W^ 



and Wi{a) = a^Wi"^^ + .... Cyclic symmetry 



implies that under a cyclic rotation P of the external legs, i i + 1, the Wi permute into 
each other according to FWi = W2, F'^Wi = W3, F^Wi = Wi, and similarly for the Wi. 
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We recall from section [21 eq. (12. 26 p . that the ratio function(s) Vi and Vi are defined 
by [360 



^NMHv(a) = -^MHv(a) [[(2) + (5)]Vi(a) - [(2) - (5)] t>i(a) + cyclic 



(6.^ 



Based on the universality of infrared divergences, and in particular the independence of 
infrared divergences on the helicity configuration, the ratio function V defined in eq. (12. 5 p 
is expected to be infrared finite, and independent of the regularization scheme used to 
compute it. More explicitly, comparing eqs. (16. ip . (16. 7p and (16. Sp . we see that 



W^{a) = M{a)Vi{a), mia) = M{a)Viia) 



1,2,3. 



(6.9) 



Expanding these relations in the coupling constant, we find, at the one- and two-loop orders, 

(6.10) 
(6.11) 
(6.12) 

(2) 

It will be a non-trivial check of our calculation that all infrared divergences cancel in . 







- M(i) , 




= w?^ 


- M(2) - M(iV/') 




= w?^ 





6.3 The one-loop ratio function 

At one loop, the MHV amplitude is given by [65] 



M(i) = -- y 

o-e5iUP5iUP25j 



F^'"((t) - ^F2"^^(a) 



+ 0(m2), (6.13) 



where = {(123456), (321654), (456123), (654321)} and (abcdef) denotes a permutation 
of the external momenta. In the NMHV case, we have 



Wi^^ = --Yl [^'"(^) + F^"''{(^)] + 0{m^) . (6.14) 

aeSi 

In writing eqs. (I6.13P and (I6.14p . we converted the corresponding expressions in dimensional 
regularization to mass regularization. The definitions of the integrals F^"^, _p2me p2nih 
in this regularization are given in appendix O 

Inserting these results into eq. (I6.10p to obtain Vi'^\ and then applying the permutation 
p2 to get V^^^ = V^^\u,v,w), we recover the expressions for l^(^) and V'^^^ in eqs. (I2.3ip 
and (I2.32p . These results are in perfect agreement with the results of an earlier com- 
putation using dimensional regularization [12], confirming the expectation that the ratio 
function should be independent of the regularization scheme. (The result of the original 
calculation [36] of the one- loop ratio function differs by a convention-dependent constant.) 



''The original definition [36] differs from one used later [42] by a (coupling-dependent) constant. We use 
the latter definition [42] because it makes the collinear behavior of the Vi simpler. Note that Vi is called 
a in ref. 
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Let us check that the colhnear and spurious conditions reviewed in section [2] are satisfied. 
They are given by eqs. fl2.38p and fl2.36p . respectively. Indeed, we have that 

hm \V^^\u,l-u,w) + V^'^\w,u,l-u)] =0, (6.15) 

and 

V^'^\u, V, 1) - V^^\l, u,v) = 0. (6.16) 
6.4 The two-loop ratio function 

There exist several representations of M*-^-* and Wi'^'^ in terms of loop integrals. Using 
generalized unitarity and dimensional regularization, representations for the loop integrand 
of M^^^ and the even part of W^^^ were found in refs. pO] and [12] , respectively. Alternative 
expressions for a four- dimensional integrand were derived using on-shell recursion relations 
in refs. [571 129]. This loop integral representation also describes the odd part W^]^^''. However, 
it will be convenient for us to choose a form in which the MHV and NMHV amplitudes are 
treated in a uniform way [TUl 112] • 

As in the one-loop case, we will assume that the loop integrals appearing in the massive 
regularization are the analogs of those appearing in dimensional regularization [12]. A 
similar assumption was made for the four-point amphtude up to four loops [751 E3], and 
for the two- loop MHV amplitudes up to six points [20]. The latter work also required 
promoting the planar four- dimensional loop integrands of ref. [57] into objects that can be 
integrated to give a finite result. We should point out that this procedure could in principle 
miss terms whose integrand vanishes as the mass vanishes, — )■ 0, but that are finite 
after integration. Although examples of such integrals have been given [751 ES], they have 
not yet proved relevant in a practical calculation. In principle, there are various ways of 
introducing mass regulators, which differ in how masses are given to different propagators, 
leading to different results after integration. We will use the mass regulator of ref. [21] . 
which provides a systematic way of introducing the masses. 

We should also comment on '/i-integrals' present in dimensional regularization, in which 
numerator factors involve explicit factors of the extra-dimensional components jl of the loop 
momentum. These integrals do not seem to have an analog in mass regularization, at least 
when one neglects terms that vanish as — )■ 0. The //-integrals arise in dimensional 
regularization due to a mismatch in dimension between the four-dimensional external po- 
larization vectors and the D-dimensional loop integration variable. It has been observed 
in explicit computations that in the quantity logM the //-integrals only contribute at 0{e) 
in dimensional regularization. At two loops, this requires a cancellation involving one- and 
two-loop /i-integrals [iQ\- Such an interference has no analog, at least through (9(m^), 
in the massive regularization, and therefore we drop the /i-integrals in the dimensionally- 
regularized integrands of refs. [TU] and [12]. 

At two loops, both the MHV amplitude and the even part of the NMHV amplitude can 
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Figure 1: Two- loop integrals J^*^ entering the two-loop six-point MHV and NMHV ampli- 
tudes. The labels i are to the upper left of each graph. Solid internal lines indicate scalar 
propagators, while numerator factors (if any) are shown below the graph. The arrow on 
the external line indicates leg number 1. The figure is from ref. |42] . 

be parametrized by |12] 

+ ici4/(^^) + ^Ci5/(^^) + Ci6/(^^) . (6.17) 

The integrals /*■'-' that enter are depicted in fig. [TJ We recall that P denotes a rotation of 
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the external momenta by one unit. The coefficients q are given by 



Cl — 


■3123 l^'^ 12645"' 234 + ■'23''56"'345 
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(6.18) 



for the MHV case [TO], and by 



Cl 

C3 
C5 

C7a 

Cg 

Cll 
Cl3 
Cl5 



"-^123-^34 -561 + -5123-5234-5345 

~-5l23-5234-5l2-545 ~ -5123-5345-523-556 +2S12S23S45S56 
-5123 (-5123-5345 ~ -512-545) 
— -512-5123-5234 

— -5 123 (-5345-5234 — -561-534) 

— -5 123 (-5234-5345 ~ -561-534) 

-5123-545-556 

~-556 -561-5123 
o2 



'123 



-561 







C2 = 2Si2-523 



for the NMHV case [12]. Here 
stood to be defined modulo 6. 



C4 

C6 
C76 

C8 
ClO 

C12 

Cl4 
C16 = 



S123-534 

-561-512-5123 

2S123-534-561 



-556-5123-5345 

~ -5 1 23 ( -5 123 -5345 ~ -5 12 -545 ) 

2-534S123 

2S12S34S123 

(6.19) 

with all indices under- 



Then we can write 



1 
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creSiUPSiup2 5i 



(6.20) 

(6.21) 



In refs. [T0| H2] . the dimensionally- regularized version of the above formulas was used to 
study these amplitudes numerically. In particular, the dual conformal invariance of the 
remainder and ratio functions was tested. The individual integrals are rather complicated, 
especially the ones of double-pentagon type, and an analytic formula for them is not known 
yet. 

Let us discuss several strategies that might be used to simplify the calculation. 

In ref . [20] , the calculation of M^"^^ in the massive regularization was simplified by going 
from the above integral basis to a more convenient one. In particular, the complicated 
double-pentagon integrals were replaced by other double-pentagon integrals (plus simpler 
integrals) that are conceptually and practically easier to evaluate. 
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Another possibility is to exploit the fact that the ratio function is dual conformally 
invariant, although the individual integrals contributing to it are not. This fact can be 
used to simplify the expression for the ratio function, by taking limits that leave the cross- 
ratios invariant, but simplify the individual integrals. This technique turned out to be very 
useful in computing the Wilson loops dual to MHV amplitudes 



Here we use a trick that relies on the following observation. The ratio between the 
coefficients cu and C13 is exactly the same in the MHV and NMHV case — see eqs. fl6.18p 
and fl6.19p . There is still a small mismatch in those terms when comparing eqs. fl6.20p 
and (I6.2ip . due to the different permutation sums. However, this mismatch disappears if 
we choose a symmetrical kinematical configuration. We can choose, for example, 

K = {xl^2 = 1 , = , ^ = 1, 2, . . . , 6, (6.22) 

which corresponds to setting all three cross-ratios equal to u. As we will see, this kinemat- 
ical subspace is more than sufficient to fix the remaining ambiguities of the ansatz in the 
preceding section. 

For equal cross-ratios, taking into account the prefactors and different numbers of per- 
mutations in eqs. fl6.20p and fl6.2ip . we see that the sum of W^'^^ and | M^"^^ not only cancels 
the contributions from I^^'^^ and but cancels or simplifies several other coefficients as 
well. We can write 

) [K] = Si^^ - |m(2) [K] , (6.23) 

where Si^^ is defined according to eq. fl6.17p . with the new coefficients 

1-u 1 111 111 . ^ 

1, 2, — ^, -, -1, 0, - - 1, 1, 0, 0, 0—,0—\, (6.24) 



and where we have combined C7 = c^a + cyh + cyc, because the corresponding integrals 
are equal at the symmetrical point {u,u,u). Given the known analytical result for M^^\ 
we only need to evaluate the integrals for i = 1,2,3,4,5,7,8,9,10,14,16 in order to 

('2') 

ohtamWl^'[K]. We could even further simplify the latter integrals using a more convenient 
integral basis |20l [STJ |29], but this turns out not to be necessary for the present purpose. 

Taking into account eq. (16. lip , we have 

\/(2) [K] = 5(2) - -M(2) [K] - M(i) [K] 1^(1) [K] . (6.25) 
3 

Let us collect the relevant formulas here, using in particular eqs. fIC.lOP and fIC.lip . and 
letting L = logm^: 

V^'^IK] = Uog'u + hi2{l-u)-C2, (6.26) 

MW[ir] = -^L^ + ^-^log^n-^Li,(l-n), (6.27) 
M(2)[K] = (logM)(2) + 1(M«[;^])2 

= lc2L'-3C3L + lu-C2F^'\K]+nf\u,u,u) + ^iM^'\K]f, (6.28) 

F^'^K] = ^-^log'u-^U^il-u) . (6.29) 
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We wish to emphasize that all terms appearing on the right-hand side of eq. fl6.25p have 
infrared divergences in the form of powers of L = log(m^), and that those terms must 
cancel in the infrared- finite quantity V^'^\ This cancellation is a non-trivial check of our 
calculation. 

The evaluation of the loop integrals proceeds in the standard way. We give a detailed 
example in appendix |Dl We derived Mellin-Barnes representations for all integrals, and 
then used the Mathematica code MBasymptotics.m [TTJ [7H1 EO] in order to perform the 
asymptotic — )■ limit. In this way we could verify the cancellation of the infrared 
divergent terms, analytically at the L^,L^ level, and numerically at the L^,L level. The 
remaining finite L° terms are given by at most four-fold Mellin Barnes integrals, which 
gives us a convenient way of evaluating V^'^\u,u,u) numerically. 

We can do better and use MBasymptotics.m another time in order to compute analyt- 
ically the small u and large u limits of V^'^\u,u,u). Having in mind that we want to fix 
the remaining undetermined coefficients of our ansatz from section [3l we go beyond the 
logarithmic terms in the expansion and also keep power suppressed terms in u. 

To promote these asymptotic limits back to a full function, we make the analog of 
the ansatz of section |3l by reducing the result of section [5] to the case of all cross-ratios 
equal. Hence we expect V^'^\u,u,u) to be given by a linear combination of TZ^'^\u,u,u) 
and single- variable harmonic polylogarithms. 

We can then compare the asymptotic limits we computed against the corresponding 
expansions of our ansatz. In fact, when fitting a complete function against just a few 
parameters it is highly non-trivial that we do find a solution. From comparing the first 
terms in the small u and large u expansion (see appendix |D] for more details), we find 



V^'\u,u,u) = -l7^f + -^log^« + 
6 lb 



Hi - jC. 



+ 2 



HI, log u + HI, + ff2,i,i - 3 & m + -3- C« ■ (6-30) 



iff," + 

17 



log^ u-H'^\ogu+ ^Hl + ^iJ2"2 



We also performed numerical checks of this expression at intermediate values of u. By 
comparing eq. (16.301) for V'^'^\u^u^u) = V{u,u,u) with our ansatz (15. 7p for v = u and 
= u, we find that the remaining twelve parameters in the joint ansatz for V{u, v, w) and 
V^'^\u, V, w) are all fixed. (As mentioned in section[5l there is one additional constraint from 
the beyond-the-symbol spurious-pole constraint, which is compatible with this solution.) 
The values of the parameters are, 

. . . ,a9,Ci,C2} = {|, i, |, -|, -i, i, 0, |, -1, 1, -1} . (6.31) 

We present the final form of the functions V and V in the next section. 



7 The final formula for the two-loop ratio function 

Now we insert the values of the twelve parameters that were fixed in the previous section 
into our ansatz, and convert everything except il*^^) and TZ^q^ into classical polylogarithms 
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whose arguments are simple, rational functions of m, v and w. The result is 

V{u, f , w) = V^{u, V, w) + V^{w, V, u) + V'^{u, f , w), 



where 



and 

V^{u,v,w) 



(7.1) 



A \ u 
1 , fuw\ 
- 4 ( VJ 

+ ^Li2(l-t;) 

+ ^112(1 -m) 



Li4(l — u) + logM Li3(l — u) 



Li3 ( 1- ^ ) +2Li3(l-n) 



Li2(l — -u) + logM logf 
2Li2(l — m) — log^f — log^w + 4 logf logw — 12^2 
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The function fi*^^) can be evaluated as a simple one- dimensional integral over classical 
polylogarithms with rational arguments, using eqs. (14. 6 p and (I4.18P from section |H The 
function TZq is the two-loop remainder function. It can be expressed entirely in terms 
of classical polylogarithms whose arguments involve square-root functions of the cross ra- 
tios Alternatively, it can be expressed, using eq. (I4.22p . in terms of three cyclic 
permutations of VL^'^\ plus classical polylogarithms with rational arguments. It is clear 

('2\ 

from eqs. (17. 2p and (17. 3p that V{u,v,w) is real in the positive octant, given that TZq and 
are. 



For the odd part we find, using ax = as 



V{u,v,w) = l{yx + f) 



(7.4) 



This is exactly the linear combination of Vx and / (multiplied by an overall |) for which 
we derived a simple parametric integral formula in section HI 

We can give an alternative form of the final answer that involves both fi*^^^ and the 
double-pentagon integral with 'mixed' numerator, fi^^-* [57], where the latter integral is 
evaluated in appendix [El Due to the respective symmetry and antisymmetry of V and V 
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Figure 2: Representation of V in terms of the finite, dual conformal loop integral 
This integral is evaluated in appendix |El The sum corresponds to the right-hand side of 
eq. (17. 9p . which also provides the proper overall normalization. 



under exchange of their first and third arguments, eq. fl2.28p . the NMHV ratio function is 
entirely specified by V{u,v,w) + V{yu,yv,yw)- We have 



V + V 



v,w) 



(7.5) 



where T{u,v,w) is implicitly defined by eqs. fl4.22p . (17. ip and (lE.7p . Explicitly it is given 
by, 

T('u, t>, w) = T"^(u, V, w) + T^{w^ V, u) + T^(u, f , w) , (7.6) 
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We see that T is given by sums of products of logarithms and polylogarithms with arguments 
which are rational combinations of u,v,w. In other words, the most complicated piece of 
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V +V is captured by the two double-pentagon integrals on the right-hand side of equation 

Moreover, the only term containing parity-odd pieces on the right-hand side of fl7.5p is 
We can easily project out the parity-even piece by taking a linear combination of this 
integral minus the same integral rotated by three steps in the twistor variables. This means 
that we have an extremely simple representation of the parity-odd function V in terms of 
finite, dual conformal loop integrals (see fig. E]), 



n^''\y^,yu,yv)-^^^\l/y^^,l/yu,l/yo)\ • (7.9) 



The same double-pentagon integral with mixed numerator appears in the representation of 
the NMHV loop integrand that was given in Table 1 of ref . [57] . The latter integral contains 
both an even and an odd part, although it is not immediately obvious how to separate the 
two. For example, although the penta-box integrals appearing in the representation of 
ref. |i57j of that amplitude contain odd parts, it can be shown that the latter are only 
(9(m^) when the integrals are evaluated using a massive regulator [21]; see ref. 



We can perform a numerical check of our result for the (parity-even) x (parity-even) 
part. Using the values obtained for fi*^^) in eqs. fl4.19p . fl4.20p and (14.2 ip . we find that 

[V + n'i'^]{f,^,^) = 14.428955293631618492, (7.10) 
[V + n'i^]{^,§,f) = 12.613874875030471932, (7.11) 
[V + n'i'^]m,f,^) = 11.705797993389994692, (7.12) 

in agreement with the values given in Table I of ref. |;42j, to the numerical accuracy given 

(2) 

there. For reference, we also give the value of TZq , which is the same for all three points 
due to its symmetry, 

n^^\f,^,§) = -3.655432869447587985. (7.13) 



We also give the numerical values of the parity-odd function at these three points. Here 
we have to specify the y values, or equivalently the branch of the square root of A that 
we consider. At the three points, A is negative, A = —1.1049134948096. We take the 
positive imaginary branch of the square root, a/A = 1.0511486549530 i in defining the y 
values through eq. fl2.17p . We then evaluate eqs. (14.330 and (I7.4p to obtain, 

0.090538030916462016642, (7.14) 
0.12117656112226985895 2, (7.15) 
0.0306385302058078423H . (7.16) 

Note that these three values sum to zero. 

In fact, although it is not apparent from the integral form (14.330 . for general kinematics 
the function V obeys 

V{yu,yv,yw) + V{y^,y^,yu) + V{yyj,yu,yv) = 0. (7.17) 
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This relation is a consequence of our ansatz and the symmetry condition (12. 28 p . Given this 
symmetry condition, eq. f l7.17p means that the totally antisymmetric part of V vanishes. 
Even if there had existed functions within our ansatz with a totally antisymmetric part, 
we could have removed them simply by noting that they never contribute to the ratio 
function due to the condition fl^Tij) . 

In an auxiliary plain text file accompanying this article, we provide the degree-four 
symbols for the functions V, V, T and Y. In these files, a term a®6®c®(i is 

written as SB (a, b, c, d). 

8 Conclusions and outlook 

In this paper we have obtained the full analytic result for the two-loop ratio function 
in planar JV = 4 super Yang-Mills theory. Our method assumed the existence of two 
pure functions, V and V, characterizing the ratio function, and was based on making an 
ansatz for the letters entering their symbols. We then further restricted the ansatz by 
imposing physical constraints, such as the behaviour in collinear and spurious regimes, 
and constraints coming from the operator product expansion of Wilson loops, leaving only 
a small number of undetermined parameters. The remaining parameters were fixed by an 
analytic computation of the loop integrals that contribute to the ratio function in particular 
kinematical regions. 

We analysed the constraints in the collinear and spurious pole limits. It is interesting 
that the spurious pole constraint involves both the (parity-even) x (parity-even) and the 
(parity-odd) x (parity-odd) part of the ratio function. We found that, within our ansatz, 
the (parity-odd) x (parity-odd) part is uniquely fixed by the (parity-even) x (parity-even) 
part. In particular, it is necessarily non-zero. 

We were able to express the ratio function in terms of sums of products of classical 
polylogarithms of rational arguments, plus two relatively simple new functions. The first is 
the parity-even double-pentagon integral fl^'^\ The second is a new function V describing 
the parity-odd sector, but it is also related to the parity-odd part of a second double- 
pentagon integral, Cl^'^^ Neither of these two additional functions can be expressed in 
terms of classical polylogarithms; however, we have provided simple parametric integral 
formulas for them, based on the differential equations that the integrals obey. We have 
checked our result for the (parity-even) x (parity-even) part of the ratio function by an 
analytic two-loop computation (in a special kinematical regime) performed in the present 
paper, as well as against numerical values in the literature. 

Let us comment on the class of functions that can appear within our ansatz. We 
considered symbols that are built from the set of nine letters {u,v,w, 1 — m, 1 — f , 1 — 
w, Uu, Vv, Vw}, with the physical constraint that the first entry should be drawn from the set 
{u,v,w} only, to exclude non-physical branch cuts. At degrees 1, 2, 3 and 4 there are 3, 
9, 25 and 69 integrable parity-even symbols of this kind. At degree 3 and 4 there are also 
1 and 6 parity-odd integrable symbols, respectively. As a byproduct of our analysis, we 
have a complete basis of functions corresponding to the parity-even symbols through degree 
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four, without imposing any symmetries or coUinear or spurious pole constraints. Three of 
the degree-four functions are given by in its three orientations, while the remaining 
functions are simple sums of products of single-variable harmonic polylogarithms, such as 
-f^o,i, 0,1(1 — u). The labels and the argument are chosen such that only a physical branch 
cut starting from m = is present. In general the labels can be any combination of zeros 
and ones, provided that the last label is 1. The unique parity-odd function at degree three 
is just the (rescaled) six- dimensional hexagon integral $6? whose relevance for scattering 
amplitudes in A/" = 4 super- Yang-Mills theory was suggested earlier [59]. The six parity-odd 
functions at degree four are the three functions ^^logu, ^glogv and ^Qlogw; two more 
functions are given by V in two orientations (which is also described by the parity-odd part 
of the two- loop mixed hexagon fi*-^-*); and there is one further function. 

Beyond two loops {i.e. for symbols of degree higher than four) new functions can ap- 
pear, as in the three-loop MHV remainder function [22j. It would be very interesting to 
find representations for them, analogous to the simple parametric integral representations 
obtained in this paper. 

The ansatz we made for the symbol was motivated by explicit results for loop amplitudes 
fl6\ [22] and loop integrals [581 El] . Another motivation comes from thinking in terms of 
twistor-space variables. Our ansatz implies that the letters of the symbol factorise into 
four-brackets of momentum twistors. This seems natural, because for six points (and hence 
six twistors describing the scattering data) intersections of lines and planes in twistor space 
always factorise into twistor four-brackets. At any rate, it would be very interesting if one 
could prove or disprove our ansatz for the six-point remainder function and ratio function 
at an arbitrary number of loops. If the ansatz is valid to all loop orders for six-point 
amplitudes in A/" = 4 super Yang-Mills theory, then it is an extremely powerful constraint 
on the S matrix of that theory. 
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A Pure functions and symbols 



We define a pure function of degree (or weight) k recursively, by demanding that its differ- 
ential satisfies 

df^'^ = J2f!^''~'^dhg<j>r. (A.l) 

r 

The sum over r is finite and (pr are algebraic functions. This recursive definition is for all 
positive k; the only degree zero pure functions are constants. The definition (]A.1|) includes 
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logarithms and classical polylogarithms, as well as other iterated integrals, such as harmonic 
polylogarithms of one [70] or more [80| [8T| [72| [82] variables. 

The symbol [T71 [TSl US] of a pure function / is defined recursively with respect 

to eq. (jXi]), 

5(/W) = 5^5(/('^-^))®0,. (A.2) 

r 

If we continue this process until we reach degree 0, we find that S{f^''^) is an element of 
the /c-fold tensor product of the space of algebraic functions, 

'5(/^')) = $^0a,®...®0a,, (A.3) 

a 

where a = {ai, . . . , ak}- The symbol of a function loses information about which logarith- 
mic branch the function is on. It also does not detect functions that are transcendental 
constants multiplied by pure functions of lower degree; such functions have zero symbol. 
The symbol therefore corresponds to an equivalence class of functions that differ in these as- 
pects. Nevertheless, the symbol is extremely useful, because complicated identities between 
transcendental functions defined by iterated integrals become simple algebraic identities. 

If a symbol can be expressed as a sum of terms, and all entries in each term belong 
to a given set of variables, then we say that the symbol can be factorised in terms of that 
set of variables. In this paper we have assumed that the pure functions associated with 
the NMHV six-point ratio function can be factorised in terms of the set (13.11) . From the 
definition of the symbol, a term containing an entry which is a product can be split into 
the sum of two terms, according to 

. . . ® 0102 8)... = ...® 01®... + ...®02®.... (A.4) 

Performing this factorisation is usually necessary to identify all algebraic relations between 
terms. It is often necessary to perform the step again after taking a kinematic limit, because 
the algebraic relations in the limit are different than for generic kinematics. 

The elements of the symbol are not all independent, but are related by the integrability 
condition d^/^'^^ = for any function f^''\ The integrability relations can be described 
simply: Pick two adjacent slots in the symbol 0q,- (g) (poi+i and replace the corresponding 
elements by the wedge product d log 0q,- Ad log (pat+i in every term. The resulting expression 
must vanish. 

The symbol also makes clear the locations of the discontinuities of the function. If 
S{f^''^) is given by eq. flA.3|) . then the degree k function /'^'^^ has a branch cut starting at 
(pai = 0. The discontinuity across this branch cut, denoted by A^^^f^''\ is also a pure 
function, of degree [k — 1). Its symbol is found by clipping the first element off the symbol 
for 

5(A^„J('=)) = X]0«,®...®0„,. (A.5) 

a 

It is instructive to check, for example, the vanishing of the double v discontinuity for the fi 
functions in eq. (13. 7p . by inspecting their symbols. Using S{H2) = — f ® (1 — f ) is enough 
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to show that /i through obey this relation. Using iS(2i7|;^ + log v H^) = —v^{l — v)^v 
and iS(if2 2) = ® (1 — f ) ® f ® (1 — f ) is enough to establish it for /s, and so on. 

In general, taking discontinuities commutes with taking derivatives, and both operations 
can be carried out at symbol level. These facts make it straightforward to verify, starting 
from eq. (13.61) . that the double v discontinuity of Vx/(2356) is annihilated by the operator 
V defined in eqs. (fSJ]) and M . 



B Details of the collinear limit 



We give here beyond-the-symbol completions of the functions Vx, fi, ■ ■ ■ , fs obeying the 
collinear limit constraint. We denote the completed functions by Fx = Vx + Vx or Fi = fi + 
fi, where the fi were given already in the main text. The coUinearly-consistent completions 
of the functions Vx and /i, . . . , /y are simple to calculate. We find that we can choose 
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(B.3) 
(B.4) 
(B.5) 
(B.6) 



log('uw) logf — log-ulogw + 2(112(1 — m) + Li2(l — w) ) + 2^3 log(-uf w) , (B.7) 



41og(Mw) logf — 9 log M log w + 8(112(1 — u) + Li2(l — w)) + C,'i\og{uvw 



(B.8) 



To define /s, the limit w — )■ of ^'•^^(w,^, 1 — u) is required. Analyzing the symbol of 
QP'\w, u,l — u), one expects the following behavior as w — )■ 0, 



lim f2*-^-*(w, u, 1 
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\og^wq2{u) + log wqs{u) + qi{u) + 0{w) . 



(B.9) 



From the symbol of VL^'^'^ we can determine the symbol of the qiiu). Therefore, the only 
ambiguities to be fixed are beyond-the-symbol terms in the gj, for which we can make an 
ansatz. Then, we fix the latter by comparing against the asymptotic w — )■ limit of a 
Mellin-Barnes representation of We find. 



q2{u) = ^log^M+ ^112(1 -m) 



(B.IO) 
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g3(^^) = -Li2(l - u)(\ogu + log(l - u)j - log^ ulog(l - u) + (2 logw 

+ Li3(l-M)-Li3(M)+C3, (B.ll) 



1 3 
Qiiu) = - log^ u log(l -u) + - log^ u log^(l - u) 
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-3S'2,2(m) +3C3logM+ -C4- 

Here 6*2,2 = -f^o,o,i,i('w) is the Nielsen polylogarithm. 



(B.12) 



The other hmit that is needed in eq. (12.381) can be obtained by the symmetry of Q^"^^ in 
the first two entries, 

^'•^^(1 — u, w, u) = Q^'^\w, 1 — u,u) = log^ w ^2(1 — u) + log w ^3(1 — u) + ^4(1 — u) . 

(B.13) 

Using these limits, we can determine a correction to /g such that /g + /s satisfies eq. f l2.38p . 

k 
3 
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log('uw) logf — logM logw — Li2(l — m) — Li2(l — ) + Cs log('Uf w) • (B.14) 



We found the following identity helpful, 
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We also have 
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because /g = TZ^q^ vanishes in all collinear limits. 

C One-loop integrals in massive regularization 

All integrals in our paper are given in the mostly-plus metric, so that the distances x]j are 
positive in the Euclidean region. 



The integrals appearing in the one-loop MHV and NMHV amplitudes are 
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Figure 3: One- loop box integrals appearing in MHV and NMHV amplitudes. 



where we recall that xf^j^^ = with indices defined modulo 6, such that in particular 
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0. See fig. [3l It is convenient to define the dimensionless functions, 



They are given by 



-ilm 



p2me 



^2 rim 
„2 Jl 



2me 



2 2 T2mh 
•^14"''26-' 



2 

9 m 


+ lo 


2 




2 

m 


g — 


g — 

^15 


-log' 


^14 


2 Lis ( 




^14 \ 


-2 Lis 


{■ 













2 ^15 



46 



14 



46 



(C.4) 
(C.5) 
(C.6) 



(C.7) 



and 



7i2me 



m 



m 



m 



m 



log' - log' + log^ ^ + log^ ^ + log' ^ 



X 



14 



36 



X 



13 



46 



X 



36 



13 



X 



14 



X 



+2Li2 (1-^1 +2Li2 ( 1-^ 



X 



2Li2 ( 1-^ 



+2Li2 1 - 



X 



46 



X 



36 



2 Li2 1 - 



^36 

"''13"''46 \ I /'o^™2^ 



X 



14 



^2 ^2 
Xi4X3g 



j2mh 



-2 Li, 



2 2 2 
^24X45 

2 N2^2 
I4J -^26 



X- 



X 



X 



- log' 1^- log , 



46 



X 



14 



X 



14 



^2 
X]_4 



l-|i)-2Li2(l-^)+OK). 



'14 



In the symmetric kinematics (16.221) . and neglecting the (9(m') terms, we have, 
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D Description of the two-loop computation 



Let us illustrate the analytic computation of the loop integrals by using the pentabox 
integral /^^°^ of fig. [H It is defined by 



/(lO) 



{m'^Y (a^Qj + ni?){x\^ + m?){xl^ + ■m?)xfj 

{xlj + m'^){xlj + m?){x\j + m?){x'lj + m?) ' ^ 



where the external dual coordinates are in the order X2, 0:3, X4, X5, Xg, Xi, reading counter- 
clockwise from the external momentum ki in the figure. Also, Xi is the dual coordinate 
for the box, and Xj is the one for the pentagon. We remind the reader that the specific 
mass assignment in eq. ( ID.ip . in particular the fact that the internal propagator is mass- 



less, follows from extended dual conformal symmetry. See refs. [21] and [73j for further 
explanation. 

We proceed by deriving a Mellin-Barnes (MB) representation for this integral. This is 
done by first introducing Feynman parameters in order to carry out the four- dimensional 
loop integrations. Subsequently, MB parameters are introduced to factorize the Feynman 
denominator, after which the Feynman integrals can be done trivially. Experience shows 
that it is convenient to introduce the MB parameters loop by loop [83j. Very detailed deriva- 
tions of MB representations for integrals like J*^^'^-' in eq. (ID. II) can be found in appendix A 
of ref. 



In the case of integral the numerator factor {x\j + m?) deserves a comment. We 
choose to treat the latter as an inverse propagator. In doing so, some of the formulas we need 
to use, such as the Feynman parameter formula, develop spurious divergences. In order to 
be able to still use these formulas, we work with the analytically continued integral l'^'^^\5)^ 
where the numerator factor is replaced by {x\j + rn?Y^'^^ and the integration measure is 
changed to d'^^'^^Xj. We will do our computation for 5 7^ 0, where all manipulations are 
allowed, and take the 5 — )■ limit later. The MB representation we find in this way is 
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(D.2) 



where 2^4^12 = X]j=4-^i- Here the integrations go from —zoo to ioo in the complex plane. 
The real part of the zi must be chosen such that the arguments of all F functions have 
positive real part. One finds that this is only possible for 5 7^ 0. 
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The limit 5 — )■ is very similar to the regulator limit in dimensional regularization, 
with the difference that here we expect a finite result, because the original integral was 
well-defined for 5 = 0. In order to take the limit, one first has to deform some of the Zi 
integration contours [H3]- This procedure has been implemented in the MB.m Mathematica 

code miEsiEn]. 

Having removed the auxiliary parameter S, we have a valid MB representation for /*^^°^ 
We can now perform the regulator limit — > 0. This again involves deforming the 
integration contours, such that the real part of the exponent of becomes positive, at 
which point a Taylor expansion in is possible. We neglect power-suppressed terms in 
m^, since we are only interested in the logarithmic infrared divergences and in the finite 
part. In deforming the contours, one picks up residues from poles of the F functions, which 
can produce powers of logm^. The resulting lower- dimensional integrals are treated in the 
same way. 

In fact, the leading divergent log*^ and log^ terms are obtained in this way without 
any remaining MB integrations. For example, 



/(lO) 



,-yi2 ^2 /-y»2 

"''24"''26"''35 



log^ + C(log^ 



(D.3) 



All log* terms with i > eventually cancel in the definition of the remainder function. 
We will therefore focus on the finite terms as — >■ 0. The latter are obtained as at most 
four-fold MB integrals. 



In the main text, we have considered the special kinematical regime K in eq. fl6.22p . in 
which all three cross-ratios are equal to u. It is easy to use the Mathematica codes [77 ] i75 | [75] 
in order to compute the m — or m — )■ oo limits of I^^^^ [K] analytically. For example, we 
find, in the small u limit, 



lim/(i°)[ir]|i„,o 



3 , 4 

32 " 



+ log^ u 
+ 0(log^ u 



1^1/2 + ^ + A^3/2 ^ 3 2 ^ 5/2 10 3 

12 36 2 12 3 ^ ^ 



(D.4) 



It is straightforward to obtain higher orders in these expansions, either analytically or 
numerically to high precision, but we refrain from reproducing them here to save space. 

(2) 

Computing the asymptotic expansions of all integrals contributing to in this way, 
we obtain 



5 , 4 

— log u 
32 ^ 

+ log^M 



+ log^M 
+ C(logM 



71 
16 



TT^ 7 2 19 3 653 4 995 5 



(D.5) 
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in the small u limit, and 



32 

— log^ u 



24 48 72 96 120 ^ 



^ M ^ M ^ U ^ U ^ U ^ + Oiu 

24 16 64 1440 8960 16800 ^ 



+ log^ u 

+ C(logM), (D.6) 



in the large u limit. We remark that the half- integer powers appearing in eq. ( ]D.4p have 
cancelled in the sum over all integrals contributing to SI . Higher-order terms in the 
expansions can be obtained numerically to great accuracy, but are not displayed for brevity. 

(2), 



Comparing eqs. fID.SP and (ID.6P to eq. fl6.25p . we can {u, u, u), or 



equivalently V{u,u,u), within our ansatz. In this way we arrive at eq. f l6.30p in the main 
text. 

We can further test eq. f l6.30p by using our four-fold MB representation for V{u, u, u) in 
order to compute some numerical values at intermediate values of u. For example, we find 

V{\,\,\) = -3.49796 ±10-^ (D.7) 
y(12,12,12) = 35.56433 ± 10"^ (D.8) 

using our MB representation of V{u, u, u), and 

= -3.497905588766739, (D.9) 
y(12, 12,12) = 35.564326922499499, (D.IO) 

using eq. ( I6.30p . 

We also note that eqs. flD.7p and flD.9p agree, within the error bounds, with the numerical 
value given in ref. |12], namely Vkrv{j, j, j) = -3.502 ± 0.002. 



E Computation of the mixed numerator integral Q*^' 
Differential equation for 17^^) 

We consider the double-pentagon integral with mixed numerator [5 

d^ZABd^ZcD (4612) (2346) 13) 



U) Vv: Vw ) 



{zTT^y {AB61){AB12){AB23){AB34:) 



^ {CDi561) n (345)) 

{ABCD){CD34){CD45){CD56){CD61) ' ^ ' ' 



where (CD(561) n (345)) = (C561)(D345) - {D56l){C3A5). 
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Loop integrals of this type satisfy simple second-order differential equations [58]. The 
key point is the presence of pentagon subintegrals that are also present in Following 
ref. [58], it is easy to see that the latter integral satisfies the differential equation 

1 ^,o^ (3461) 



(2346) (1234) (2346) 

where the (rescaled) one-loop hexagon integral with mixed numerator is defined as 
(4612)(2346) f cL^Zab (A513)(A5(345) n (561)) 



(E.2) 



ui Vvi Vw) 



(3461) 



iTT^ {ABQ1){AB12){AB23){AB3A){ABA5){AB56) ' 

(E.3) 



It is given explicitly by [29] 



^l''^\yu,yv,yw) = logwlogt;- ^"^^^ ^"^ logt;logw- / logwlogw. 



(E.4) 
(E.5) 



1 - VuVv 1 - VuVv 

We note that the integrals Cl^'^^ and (l^^^ are left invariant by the transformation 

Zi < — > , Z4 i — > Zq , 
which implies 

Ui — >v, Vu — >i/yv, Vv — >^/yu, Vw — >^/yw (E.6) 

We make the ansatz that (l'^'^^ is a pure function, whose symbol's entries are drawn from 
the set of nine letters {u, v,w,l — u,l — v,l — w, y^, Vw}- Within this ansatz, we find 
that eq. (\E.2\i has a unique solution obeying the symmetry condition (lE.6p . integrability, 
and the first entry condition. The solution involves parity-even as well as parity-odd terms. 

Having determined the symbol of fi*^^^ from the differential equation (lE.2p . we now 
promote it to a function. We find that we can express it a^ 

n^^\yu,yv,yw) = ^ [n'^^\v,w,u) + Q^^\w,u,v)] +Y{u,v,w) + 2V{y^,y^,yu) , (E.7) 



with 

where 

F^(M,t;,ti;) 



Y{u, V, w) = Y^{u, V, w) + Y^{v, u, w) - Y^{u, v, w) 



(E.8) 



4Li4(M) - Li4 ( 1 



u 



+ logM 



+ 2 log(l 



u 



Lis 1 



2Li3(l -u) + 3Li3 ( 1 
2 1 



u 



+ - log U U2 I 1 



- log^u-2r(u;) + 3Li4 ( 1 - - 
6 \ w 



log ( - 



2Li3 1-M +Li3 1-- 

U 



logMLi2(l - m) - - log^M 
6 



w 



+ - log' - Li2 1 - - 



U 



(E.9) 



""To avoid confusion, we emphasize that V{yv,yw,yu) differs from ¥V{yu,yviyw), where P denotes a 
cycUc shift of ah twisters by one unit. In fact, we have P {u, w, w, ?/„,?/,„, y^} — {v,w,u, l/j/i,, l/yw, ^/Uu}, 
and hence VV{yu,yv,yw) = -V{yv,yw,yu)- 
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16 1 16 

(2) ^ (2) 

n ( yu 'Jv 'Jw) p^^o limit - a ( yu Jv 'y^) 

Figure 4: The integrals fi*-^-* and —Vl^'^^ have the same soft limit ps — )■ at the integrand 
level. This property allows us to formulate the boundary condition (lE.14p . 

with r[w) defined in eq. f l4.24p . and where the beyond-the-symbol ambiguity for a function 
symmetric in u and v is given by 

Y^{u, V, w) =C2 [ci(Li2(l -u) + Li2(l - v)) + cs Li2(l - w) + C3 (log^ u + log^ v) 

+ C4 log^ W + C5 log U log V + Cq log(Mf ) log W~\ 
+ Cs [C7 log{uv) + C8 log w] +09(4- (E.IO) 

We can ask how many of the Cj can be determined by the differential equation (IE.2p . Using 
the variables from appendix [F] it is not hard to verify that the only functions appearing 
in {u,v,w) that are annihilated by the differential operator are ^4 and (slogi^w/ (uv)). 
Therefore, 7 out of the 9 coefficients Cj can be determined by plugging eq. flE.7p back into 
eq. (Q. 

Indeed, using the parametric integrals derived in the main text for and V, we can 
easily verify the differential equation flE.2p numerically. We find 

ci = 1 , C2 = -2 , C3 = 3/2 , C4 = -1 , C5 = , C6 = , C7 = 2 - cg . (E.H) 

We will fix the remaining two free parameters Cg and Cg from boundary conditions that we 
discuss presently. 

Boundary conditions for Q^^^ 

Let us discuss appropriate boundary conditions for ^l^'^K Here we can use our previous 
experience with the integral which at the integrand level differs from fi'-^-' only by the 
numerator in one of the pentagon subintegrals. In fact, the numerators of the two integrals 
are given by 

N{n^^^) = (4612)(2346)(A513)(CD(561) n (345)) , (E.12) 
N{n^^^) = (2345) (5612) (3461) (A513)(CD46) . (E.13) 

Previously it was observed that the integrands of these two integrals reduce to the integrand 
of a penta-box integral in the soft limit ps — )■ 0, or equivalently — )■ + /JZg, as shown 
in fig. m ^20j. Unfortunately, the penta-box integral is infrared divergent, so that the limit 
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is more subtle at the level of integrals. However, we can use the fact that the numerator 
N('tl^'^^)+N{il^'^^) vanishes linearly in the soft limit. Because the explicitly known penta-box 
integral [20] only has logarithmic divergences, we expect the following boundary condition 
to hold, 

lim + r](2) ) ^2T, 1 - r) = . (E. 14) 

Here we have parametrized the soft limit for the cross-ratios u, f , w by r — )■ We have 
verified equation flE.14p at the symbol level. In the following we will assume it holds also 
at the level of functions. 

A related observation is that $7^^^ vanishes in cyclically related soft limits, 

lim fi(2)(l-r,eir,6T) =0. (E.15) 

This vanishing can in fact be understood as a property of the pentagon sub-integral. Since 
n^^) contains the same sub-integral as Q^^^ , we expect the same boundary condition to hold, 
i.e. 

lim ^](2)(l-r,elr,67-) =0. (E.16) 



We find that imposing the two boundary conditions (IE.14P and (IE.16P fixes all but one 
of the beyond-the-symbol ambiguities in eq. flE.lOp . 

ci = 1 - C9/5 , C2 = -C9/5 - 2 , C3 = 3/2 , C4 = -1 , C5 = , ce = , Cy = 2 , Cg = . 

(E.17) 

Comparing to eq. flE.llI) . we see that the two solutions are compatible with each other, 
which is a non-trivial cross check. Moreover, taken together they uniquely fix all the 
beyond-the-symbol parameters, and we have finally, 

ci = 1 , C2 = -2 , C3 = 3/2 , C4 = -1 , C5 = , C6 = , C7 = 2 , cs = , Cg = . (E.18) 



F Useful variables 

In this paper, we found it useful to work with several sets of variables. We can express the 
letters appearing in our symbols in terms of four-brackets of twistors, 

(6123)(3456) (1234)(4561) (2345)(5612) 

(6134) (2356) ' (1245) (3461) ' (2356) (4512) ' ^ ' ^ 

^ (1356) (2346) ^ _ (2461)(3451) ^ (3512) (4562) 

(1346)(2356) ' (2451)(3461) ' (3562)(4512) ' ^ ' ' 

_ (2361)(2456)(3451) _ (3462)(3512)(4561) _ (1246)(1356)(2345) 

^" ~ (2351) (2461) (3456) ' ~ (3412) (3561) (4562) ' ^'^ ~ (1256) (1345) (2346) ' 

(F.3) 



^There is a slight abuse of notation here since, strictly speaking, fj^^^ should be thought of as a function 
of the y variables. However, in the soft limit, its parity-odd piece vanishes, justifying the use of the u 
variables. 
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Since the twistors are redundant, it can sometimes be useful to have a particular parametri- 
zation for them, e.g. 



=(1, 1, 7, 1) , Z2 = (1, 0, 0, 0) , Z3 = (0, 1, 0, 0) , 

=(0,0, 1,0), ^5 = (0,0, 0,1), Z6 = (l,a,l,/3), (F.4) 

with 

a = - , P = — , J = - . F.5 

1 - VvVw i--yw 1 - VuVw 

Although the y variables are constructed using square roots of the original cross ratios u, 
V and w, the cross ratios themselves are rational combinations of the variables and 
y^. The explicit relations are, 

yuil - yv)il - yw) y^{l - y^){l - yu) |/^(1 - y„)(l - y^) 



{^-ywyu){^-yuyv) (1 - y„t/^)(i - i/^y^) (1 - yt,y«,)(i - i/«,i/«) 

(F.6) 

{l-yu)il-yuyvyw) ^ {I - y^){l - yuyvyw) .-^ „^ 

(1 - ywyu)[l - yuyv) (l - yuyv)[l - yvyw) 

1 {I - yy,){l - yuyvyw) nr (1 - y„)(l - l/^,)(l - y«,)(l - ynyt,y^„) . . 

^~^ = 71 \^^ V' = wTj X , (F.8) 

(1 - y^y^)[l - y^yu) (1 - 2/«i/«)(l - yvyw)[^ - ywyu) 



where we have picked a particular branch of v A. 
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